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PREFACE 


In recent years, the trend of airborne radars has been toward shorter 
and shorter wavelengths where the atmosphere becomes increasingly more 
opaque and weather clutter becomes increasingly more frequent and intense. 
Since either or both of these factors may seriously limit the use of such 
shart wavelengths, an assessment of their magnitudes is of considerable 
importance. A comprehensive evaluation of this sort was made by Ryde 
(1946) whose work remains as the one outstanding authoritative reference 
in this field. 

Ryde concerned himself with the extremely broad problen related to the 
detection of targets in and through clouds and precipitation. His interest 
included the wavelength region from 1 - 10 cm, with emphasis on the upper 
end of this band. Among other things, he considered the basic scattering 
and absorption cross-sections of the various atmospheric constituents and 
applied these to the computation of the magnitude of the echo intensity 
and attenuation due to each canponent individually. Although Ryde's work 
is applicable almost in its entirety today, the results of more recent 
research should be added, and minor modifications made. In addition, 
detailed consideration should be given to the wavelengths shorter than 
3 em, the region which is just now being fully exploited, Finally, it is 
of some practical interest to evaluate the cumulative effects on the 
atmosphere's transparency of the various atmospheric constituents in the 
amounts and combinations which are likely to be found under typical or 
extreme weather conditions, These are, in essence, the objectives of this 


report, 


In Part I, Marshall, East, and Gunn have made a fairly comprehensive 
survey of the state of knowledge of scattering and attenuation by the 
atmosphere. This includes some previously unpublished and important work 
by Haddock (1948). It contains in direct or indirect form the answers to 
most of the solvable problems which have been, or may be raised in this 
general field. However, an all-inclusive study would require practically 
a complete repetition of Ryde's excellent 1946 paper. Because of the 
comprehensive nature of the latter paper, it is recommended that this 
report be used in conjunction with it. Wherever the present work supersedes 
Ryde's paper, by virtue of the use of more recently determined constants 
' or theoretical information, it will be spect riest ie mentioned. 

Perhaps a statistical study as that made by Bussey (1950) would be 
of some practical utility to those interested in the average effects of 
the atmosphere on microwave transmission. In that study, Bussey computed 
the frequency with which various attenuation coefficients could be expected 
to be exceeded in the vicinity of Washington, D.C, The results depend, of 
course, on the local rainfall and water vapor statistics although they may 
be extrapolated, with caution, to areas of similar rainfall regimes, The 
information is restricted, furthermore, to the lower layers of the atmosphere, 

The present and potential uses of microwave techniques by the Air Force 
generally involve transmission through the vertical dimension of the atmos-— 
phere. Since the vertical distribution of the several attenuating elements 
(rain, snow, clouds, oxygen, and water vapor) is generally very complex, and 
since few statistics are available on the frequencies with which they occur 


in various combinations, no satisfactory estimate can presently be made of 


a3 


the amount of ee which losses of certain magnitude can be expected. 
Instead, in Part II of this report an attempt has been made to select some 
typical idealized cloud and precipitation conditions and to compute the 
resultant opacity of the atmosphere in each case, The reader may then select 
one or more of these models and corresponding losses as typical of the 

average or extreme conditions under which his radar may be required to operate. 
Such an approach to the problem may well result in a better estimate of 
typical losses than a statistical survey which is based on the necessarily 
crude use of scanty data, 

The weather models of Fart II are based primarily on physical considera- 
tions, and are therefore largely independent of geographical location. The 
frequency of occurrence of each type of situation may vary widely from place 
to place, but no attempt has been made to treat such climatological problems. 
However, the cloud densities and rates of precipitation have been described, 
by adjective ratings of "high, low, or moderate" corresponding roughly with 


conditions generally observed in temperate United States. 


An Appendix, a report entitled "Calculated Sensitivity of Airborne 
Weather Radars", became available during the publication process. Acknowledg- 
ment is made to J. S. Marshall and Walter Hitschfeld of McGill University for 
this paper. The original work was performed under joint sponsorship of 
Aeronautical Radio, Inc., and the Air Transport Association of America. 
Aeronautical Radio, Inc., has kindly granted permission for reproduction of 
this paper as an appendix to this report. 

The authors are indebted to Mr, Ralph J. Donaldson for valuable criticism 
and editing. 
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PART I 
THE MICROWAVE PROPERTIES OF PRECIPITATION PARTICLES 
by 


Je Se Marshall 
T. We. Re East 


Ke: Le Oe Gunn 


McGill University 


Montreal, Canada 


1. INTRODUCTION 


In the comparatively new subject of radar-meteorology, the theoretica]. 
paper by Ryde (1946) stands out. In it he predicted the back-scatter and 
attenuation by meteorologica]. phenomena, and, in almost every respect, the 
results obtained in experiments have confirmed the predictions. The present 
contribution will review the theoretical problem in the light of developments 
in the years that have elapsed since Ryde's work was done, and with reference . 
to the experimental results. 

Meteorological radar equipment provides information about the structure 
or pattern and about the quantity of precipitation. This paper will: not be 
concermed with the structural information obtained, but will consider the 
scattering from a small region in which the properties of the precipitation 
can be taken to be uniform, and the attenuation taken along a path of 
propagation as a whole. | 

Plane electromagnetic waves travelling through air containing precipi- 
tation are scattered and absorbed by the particles of ice, snow or water. 

Ice particles scatter some of the radiation in all directions, but in most 
cases of interest do this conservatively, without converting any into thermal 
energy, since ice is very nearly a perfect dielectric. Water, with its larger 
dielectric constant, scatters more strongly than ice; in addition, it has 
considerable dielectric loss and a large part of the attenuation by rain is 
due to thermal dissipation. 

Three quantities with the dimensions of area are derived for a particle 
in the path of a plane travelling wave. The back-scatter cross-section Tv is 


defined in such a way that & multiplied by the incident intensity would be 


ie 


the total power radiated by an isotropic source which radiates the same power 
in the backward direction as the scatterer. The "scattering cross-section" 
Q, is such that Qs? multiplied by the incident intensity, is the power 
scattered by the particle. The "total ahsorption cross-section" Q; is 

such that Q, times the incident intensity is the total power taken from 
the incident wave. 

The scattering properties of precipitation have made possible extensive 
use of microwave radar sets for meteorological purposes. The radar signal 
from precipitation is, however, not steady like that from a point target. 

The received power at any instant is made up of the resultant of the signals 
from a very large number of particles, and depends on their exact arrangement 
in space, which is being continually changed by the turbulent motion of the 
air. Marshall and Hitschfeld (1952) have shown that an instantaneous obser- 
vation of received power gives very little information about tne precivitation, 
and the signal at any given range has to be averaged until a large muimber 

of independent returns have been received. 

This can be arranged by using a long-persistence display or by photo- 
graphic recording so that the returns from many transmitter pulses are 
combined. 


The equation for the average received power Ee is 


BE Na ss 
Poe Ar Gi 
: Snr 1 


where P, is the transmitted power, 
h is the length of the pulse in space, 
r the range from which the signal is being received, 


A, the equivalent antenna aperture, 


€ 


210 the total back-scatter cross-section of the particles in unit volume. 
1 


o2u 


The received nower is proportional to h_ so that the sensitivity of 
the system can be increased bv lengthening the transmitter pulse. (If this 
change is accomranied by a reduction in receiver bandwidth, a further improve- 
ment in signal-to-noise ritio is obtained.) However, the range resolution is 
correspondingly reduced. Increasing A,, on the other hand, improves both 
sensitivity and angular resolution,-so that the antenna is nade as large as 
practicable. 

The equivalent averture A. is less than the actual area of the antenné 
for two reasons. Firstly, not all the power coming from the primary feed 
falls inside the parabolic reflector, neither is the distribution of power 
unifori across the aperture (see e.g. Silver (1949) section 6.4). In normal 
radar systems the effective aperture A! for a point target is about 65% of 
the actual aperture. The second reason which apvlies aries distrivuted 
targets affects werther radars. Some of the radiated ynower goes into side- 
lobes and the eiges of the main »beain, and thoush all of it strikes precivitation, 
scme of the returned rower comes from directions to which t*e« antenna is less 
sensitive for reception. This means that the aversse receiving apertur2 A, 
for all parts of the target is less then A!', the value for the central 
portion. 

The fact. that the target normally fills the nein at any ranzve escounts 
for the inverse square law for reccived rower versus ranze, in place of the 
usual fourth pnower law for 2 point tarzet such as an aircraft. Though 
equation (1.1) does rot mentior wavelength explicitly, @ increases ranidly 


as the wavelensth cecreases, as discused in the next section. 


=3- 


It will be shown in section 3.2 that for a wavelength 2X long 


comvared with the diameter D of the particles, the approximation 


56 | 
T = NE: \x|* peer). 


can be used, where }K]© depends on the dielectric constant of the material 

and varies only slowly with wavelength (chapter 2). From (1.1) and (1.2) we 

gee Boe TPe teh gt 16 uy? (1.3) 
ar= a4 Hi 

So fer we have neslected the attenuating effect of any precipitation 

through which the signal may have to pass. 


The intensity of plane parallel radiation in traversing distance x 


through uniform precipitation is reduced from I, to I, where 


I = I, exp {-« 4 2. | ‘ aS 


The summation is over all the particles in unit volume. In practical units, 


the attenuation per kilometre of path length is 


attenuation (db/km) = 0.4343 oP as (1.5) 
x 


where the summation is over one cubic metre and Q is now in em“. For 


the case of wavelength long compared to the drop diameters (section 3.2), 
we can use an approximation for Q., equation (3.7), so that equation (1.5) 


becomes 
attenuation (db/km) = 0.4343 um ED? Im (-K). C156) 
aN x : 


For radar, the radiation travels there and back, so the numerical constant 


should be doubled to obtain decibels per kilometre of range. 


oh 


The complete ecuation for the average received power Pel is, then 


= -0.8686 x £Q 
10 1 . 


rU 
i 


1908088 x 2 ela 


where the Rayleigh anproximation is applicable. 
If the attenuation of the precivitation is not uniform along the path, 


the factor’ ‘x rb ore must be replaced by fe Q,.dx over the path. 
iF ] 


2. RGFRACTTVE INDEX 


To calculate the scattering and attenuation of electromagnetic waves 
by particles in the atmosvhere, one must know the values of the complex 
dielectric constant € for the temperatur? and wivelenvtns under consideration. 


The square root of the dielectric constant m may be written 
Mkt — a 5s 
where n is the refractive index and K the absorption coefficient of the 
mater‘al. Jn teneral. the scatterint and attermiatior are complicated 
functions of m (chapter 3), but in the special case of wavelength long com- 


pared to diameter, 


ss 
back=-scatter of |K|” assy, 
attenuation of In (-K) b1.6) 
where K = meee : 
m* + 2 


For a given substance m is a function of wavelength and temperature. 
Saxton (1946) measured values of n and K, for water at a few wavelengths 
and calculated other values in the centimetric band, for temperatures from 
0° to 40°C. Recently Lane and Saxton (1952) have reported additional 
measurements at a lower wavelength (6.2 mm) and some more accurately deter- 
mined values at 1.24 and 3.2 cm. As well, they have determined values of k 
for supercooled water down to -8°C at three wavelengths (see section 4.3). 

In table 2.1, the variation of the refractive index with temperature 
and wavelength is shown, and the quantities Im(-K) and K @ have been 
calculated. Values for n and K are taken from Kerr (1951) for 10 cm wave- 
length, and directly from Lane and Saxton (1952) for the other wavelengths. 

Measurements of the quantities n and k, for ice show that in the 
centimeter band both n and K, are practically independent of A. Dunsmuir 
and Lamb (1945) and Lamb (1946) found k, to be a very small quantity, 
ranging from 1.2 x 1073 at OC, through 1.9 x 1074 at -200 to 1.3 x 1074 at 
-4,0C. More recent measurements by Cumming (1952) give values for k roughly 
twice those of Lamb; Cumming's value at -12C agrees with one obtained at 
that temperature by the MIT Dielectric Laboratory (Cumming 1952). Both 
Cumming and Lamb and Turney (1949) agree on a value of n = 1.78 for ice, 
constant over the temperature range 0 to -20C. The refractive index data 
for ice is summarized in Table 2.2. | 

Snow may be considered as a mixture of two dielectrics, ice and air. 
According to Debye (1929), one can calculate the ratio 7 for a mixture of 


two dielectrics by adding the ; values for the two substances in proportion 


apa 


Table 2.1 


Variation of Refractive Index of WATER with Temperature and Wavelength 


roc / A cm 10 3.21 1.24 0.62 
20 8.88 8.14 6615 Lyell 
10 9.02 7.80 5 6h5 3.94 
i ) 8.99 Tolk 4075 3045 
<8 6.48 4.15 3.10 
20 0.63 2.00 2.86 2059 
k, 10 0.90 ells 2.90 2.37 
0 1.47 2.89 Patt, 2.04 
2.55 1377 
20 928 ©9275 ©9193 «8926 
‘xl? 10 9313 9282 29152 .8726 
0 ~9340 «9300 ~9055 28312 
<8 ~8902 7921 
20 00474 .01883 20471 0915 
a) 10 .00688 0247 .0615 1142 
.01102 .0335 .0807 SE ANI 
-8 .1036 .1713 


to the mass of each. Neglecting the contribution of ; for air, then ; 
for ice of any density is practically constant. Cumming (1952) made 
measurements of the dielectric constant of snow of various densities. Using 
his data one finds that 5 varies more or less uniformly from 0.46 for ice 


of density 0.917 gm m™3 to 0.50 for density 0.22 gmom73, 


Table 2.2 


Variation of Refractive Index of ICE with Temperature 
(based on data from Cumming, 1952) 


BEE in 


7.9 x 1074 


5.5 x 1074 


0.197 At all temperatures. This is 
for ice of unit density, the 
value to be used when D is 
diameter of melted ice particle, 
in equation (1.3) Marshall and 


Gunn (1952) . 


906 x 1074 
3.2 x 1074 


2.2 x 1074 


3 SPHERICAL PARTICLES 


3.1 Mie's Theory 


A complete theory for spherical particles of any material in a non- 
absorbing medium was developed by Gustav Mie (1908). Mie's work has been 
restated by Stratton (1941), Goldstein (1946) and by Kerr (1951). Using the 
notation of these later workers, the cross-sections of a spherical particle 


are 


REID: ran A 2 
i iF Z (a) (one alata 7 b,)| (3.1) 
a eo 
a Be S(t 2) Clenl® + Veal) 0.2) 
n2 3 
and, = XX (-Re) 2 (2n + 1) (a, + b,) (3.3) 


where A # wavelength in the mediun, 
a, = coefficient of the nth magnetic mode, 


by = coefficient of the nth electric mode. 


The coefficients a, and b, are given in Appendix II in terms of 


Spherical Bessel Functions of order n. The arguments of the functions are % 
and m_ and contain the properties of the sphere and the wavelength. The 
quantity = m, where D is the diameter of the sphere. The quantity m 
is the complex refractive index of the material relative to the medium, and 
m=n- jk, where n is the refractive index, and K, the absorptive index. 
Equations (3.2) and (3.3) have been computed by Lowan (1949) for water 


spheres at 18°C for various wavelengths between 3 mm and 10 cm. 


~G~ 


3.2 Small Particles (Rayleigh Approximation) 
By expanding the a's and b's in ascending order of KA for n=l, 2..., 


formulae for Q, and Q can be obtained which apply to small values of a. 


Neglecting terms of higher than sixth power of &, only ay» b, and b, are 
significant, and they are 
a=: aL (m@-1) om? (3.4a) 
45 
pends Si me ; 
b= - A madly? (+3 2-2 «x? - 21 a «>) (3.4) 
3 m242 > m2+2 3 m@+2 
2. 
by = oy Ek 7. (3ehe) 


15 on203 


For diameters small compared with the wavelength, we can obtain approxi- 
mate formulae for QJ and Q, by putting &<1 into these equations. Then, 
only the leading term of b,, the coefficient of the first electric mode, need 
be considered, and the contribution to W7 and Q, of all higher order terms 
in the a's and b's can be neglected. The resulting formulae are know as 


the Rayleigh approximations: 


2 21012 
g =X .2«° m@=l wore [KI 6 (3.5) 
aT m*+2 rn& 
peony ee caus (3.6) 
4 8 Orr e 3 es ° e 


The total absorption cross-section Q is given by equation (3.3) in 
which the real parts of the coefficients @ and b have to be obtained. 
In this case, however, it does not always happen that the leading term in by 
contributes the most to Qe Indeed, if m is purely real, the «3 and A? 
terms of b, are imaginary and contribute nothing to a » and of the terms 


listed in equation (3.4) only the %° term of b, will appear in Q. 


-10- 


At this stage it is convenient to introduce a fourth cross-section Q,, 
the absorption cross-section, to account for the power taken from the incident 
wave which is dissipated internally as heat and not scattered; it is simply 


Q, - Q,. Putting XK1 the contribution to Q of the «° 


term in bys is 
the same as the leading term in Q,, so that this term does not appear in Q,. 
Then the leading term in Q, is given by (3.6) and the leading term in Q, 


comes from the of? term in b, and is given by 


nw; 3 ( m@-1 
Q = ,207 Im/- t= (3.77) 
eT atT =) 
and, of course, 
% =@, +Q,. (3.8) 


For a pure lossless dielectric Q. = Oso that Q =Q, and is given 
by (3.6). For a very lossy dielectric like water and particle sizes such 
that X<<1, Q5<< Q, (figure 3.2) so that Q, = Q, and is given by (3.7). 
For ice at OC and in the size range of snowflakes, at centimetre wavelengths 


Q, and Q, are comparable and both must be calculated; Q, is their sum. 


3.3 Larger Particles 

F, T. Haddock (1948) has computed for water at 18°C the ratio between 
Q, calculated from Lowan's (1949) Bureau of Standards tables, and Q, cal- 
culated from equation (3.7). This ratio Q,(Mie)/Q,(Rayleigh) is plotted in 
figure 3.1 for various wavelengths. Ryde's (1946) figure 3 is a similar 
plot, but of a, (Mie) 


_ see ee 
a geometrical cross-section 


for a wavelength a little over 1 cm. 
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FIG. 3.1 - Ratio of actual attenuation (upper) and actual back-scattering 
(lower) to that given by the Rayleigh approximation, for water 


at 18C. 
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We can interpret these curves a little further by examining the relative 
magnitudes of the terms in equation (3.1). For example, for water at 18°C, 


and A = 3.0 cm, we find that the parts of the a's and b's affecting Q 


ee Re [a,J ariepa espns 
Re [by] = (0.0121 0&7 + 0.0160 x) + 0.411 xo 
Re [b,] = - 0.000515 &? 

up to «6, 


As long as Re [b,) is much greater than Re [a,]J, only the oo? tern 
is important, but at = 0.131, Re (a,] = Re (b,J, so that Q, is about 
twice the Rayleigh value. As increases beyond this, Q tends to an 
ne law, until other terms become more important. 

The region of Q& in which the exact value of Q, given by (3.3) and the 
approximate value given by equation (3.8) [Hees by (3.6) + (3.7) | are nearly 
equal is known as the"Rayleigh region for attenuation". The diameters of 
cloud particles are well inside this region for all wavelengths down to 3 mm, 
Rain will be discussed in section 4. 

The corresponding situation for back-scatter is shown also in figure 3.1. 
Here T (MieY U(Rayleigh) is plotted against o{ for various wavelengths. We 
shall define the "Rayleigh region for back=-scatter" as the region for which 
this ratio lies between $ and 2. It will be seen that for wavelengths up to 
17 mm, the ratio stays within these limits up to about & = 1, where it falls 
rapidly. However, at longer wavelengths the ratio first rises beyond 2 so 


that the region only extends to @= 0.4 or less. 


Ete 


At 17 mm wavelength, the Rayleigh region as defined here extends to 
6 mm diameter and includes all normal raindrop sizes. At shorter wavelengths, 
the limit occurs at smaller diameters, so that at A = 3 mm, only particles 
less than 1 mm diameter are inside it. At 3 cm wavelength, the region is also 
reduced, extending only to 3 mm diameter, but at A = 10 cm it covers all 
raindrop sizes so that the approximate formula can be safely used on rain. 

The ratio Q,/Q, for water at 18°C, that is the proportion of the total 
absorption cross-section which is due to scattering, has been computed from 
Lowan's (1949) tables. This is presented in figure 3.2 for values of X up 
to 5 (i.e. covering a range of diameter up to 0.6 cm and wavelengths from 3 mm 
to 10 cm). It will be seen that the curves very roughly coincide, so that for 
diameters and wavelengths where & = 0.8 this fraction is 50%, and where 

o4 = 0.5 the fraction is 25%. There seems to be no obvious physical reason 


for this result. 
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FIG. 3.2 - Ratio of scattering to total cross-sections for water spheres. 
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3.4 Water-Coated Ice Spheres 
Aden and Kerker (1951) extended Mie's theory to concentric spheres. 
k 


Langleben and Gunn (1952) have applied this extended theory to spheres of 
ice coated with water in the region @&<1. They consider a sphere which 
changes gradually from all ice to ice coated with water, and finally to all 
water, maintaining constant mass. The calculated cross-sections for such a 
particle are plotted against the ratio of the mass of water to total mass in 
figure 3.3 for 3 cm wavelength and 1.2 mm diameter initially. As the thickness 
of the water coating increases from zero, the quantities 7, Q, and qe 
all increase rapidly from their values for ice; when the sphere is still 70% 
ice (so that the water thickness is less than one-tenth of the total radius) 
GV and Q, are already nearly equal to those for the all-water sphere. 
Earlier still, Q reaches a maximum value, about twice as great as that for 
water. No simple explanation can be offered as yet for this extremely rapid 
increase of the attenuation to twice the all-water value. For other wave- 
lengths and diameters in the Rayleigh region, the curves for Qs: voand Q 
are very similar to those of figure 3.3. For larger values of oc, the 
scattering cross-sections do not rise so early, and Qb does not show such 


an extreme behaviour. 


k 

Langleben, M. P. and K. L. S. Gunn, 1952: Scattering and absorption of 
microwaves by a melting ice sphere. Stormy Weather Research Group, McGill 
University, Scientific Report MW-5, under USAF Contract AF 19(122)-217. 
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FIG. 3.3 - Relative cross-sections of water-coated ice spheres. 


4, RAIN 


4.1 Drop Size-Distributions 


Where Rayleigh theory is applicable, radar can be used to measure the 
quantity Z& po, A more significant quantity meteorologically is the rainfall 
rate, R sie p z Yn D3 (where Vp is the terminal velocity of a drop of 
diameter D). To relate the two quantities the drop size-distribution must be 
observed. Distributions for rain at the ground have been determined in various 
parts of the world, and the distribution law, the water content M, the rain- 
fall rate R, and the radar quantity Z can all be related empirically. 

(For example, see Wexler (1948), Marshall and Palmer (1948), Hood (1950) etc.) 
Various workers have calculated simultaneous values of Z and R from the 
distribution of drop sizes obtained on a horizontal surface (e.g. Marshall, 
Langille and Palmer, 1947). Then from a number of such simultaneous values 
plotted on log-log paper a locus of best fit for the data of the form 

Z = aR? can be computed. There have been many such Z/R relations quoted 

in the literature in the past ten years, which give a wide range of Z values 
for a given R. Out of these, the authors have found eight relations which 
are based on measurements using raw drop size-distribution data. These, along 
with pertinent information, are included in Table 4.1. 

The number of points used to make up the different loci varies from 
worker to worker, as does the range of rainfall rates covered, The great 
majority of the points in each case lie between 5 and 50 mm hr7l, 


A locus 
Z = 243 pl.57 + a factor 1,22 (4.1) 


iat 


TABLE 4.1 


Zz Author Location How locus determined 


208 R153) Wexler (1948) Hawaii least squares fit (49 points) 
on Anderson's 


(1947) data 


& 
269 R1e55 Mt. Washington Cambridge, Mass least squares fit (63 points) 
observation (1951) 


214, R1le58 Wexler (1948) on least squares fit (98 points) 
Laws and Parsons! 
data 

200 RleS = Revised from Ottawa best fit by eye (with a check 
Marshall and to centre locus amongst data 
Palmer (1948) at that slope) (250 points) 

295 RLe61 Hood (1950) Ottawa best fit by eye (270 points) 

505 RL-44 Best (1950) Shoeburyness 

257 Rle4d Best (1950) Ynyslas 

1436 RLe°4 Best (1950) East Hill 


will define a band taking in the first five relations. If we include Best's 
data, the locus becomes (Z = 364 pie) + a factor 1.75). As well, one mst 
recognize that there is a limit on the accuracy with which Z or R can be 
calculated when samples are taken from relatively small volumes of space. The 
natural variations due to the methods of sampling usually employed may lead 

to an uncertainty in Z of as high as = 50%. Thus, if one uses a relation 
such as (4.1), it is probably good to within a factor 2 over a range of 
moderate rainfalls, say 0.5 to 70 mm/hr. Until more is known about the 
relationship of the radar quantity Z to a meteorological quantity, such as 


R, this uncertainty will have to be accepted. 


* Mount Washington Observatory, 1951: Contribution to the theory of the 


constitution of clouds. Quarterly Progress Report of the Mount Washington 
Observatory staff under USAF Contract No. AF 19(122)-399, Oct.23. 
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4.2 Scattering by Rain 

The curves in figure 3.1 showed the behaviour of the exact back- 
scattering cross-sections for single drops. Haddock, in plate 6 of his 
unpublished report, has evaluated the exact cross-sections for rain of various 
intensities using the drop size-distribution results of Laws and Parsons (1943). 
Haddock's data has been replotted in figure 4.1 which shows &@ against 
rainfall for two wavelengths, 3 cm and 0.9 cm. For haan? the Rayleigh 
approximate back-scattering cross-sections are shown for these same wave- 
lengths. They are derived from summing equation 3.5 over unit volume and 
“ingNtanhare ape epee oe 

Haddock Re not give a plot of the exact ZW for wavelengths greater 
than 3 cm, but the 10 cm curve will be Meters the Rayleigh one at all 
but the very heaviest rainfalls. At 0.9 cm, the exact back-scattering cross- 
section is larger than the Rayleigh for low rainfalls; it is the same at a 
rainfall of about 10 mm/hr and then becomes less quite quickly at higher 
rainfalls, being some 7.5 db below the Rayleigh value at 10 mn/hr. This 
rapid divergence of the two functions at 0.9 cm for heavy rainfalls is due 
to the rapid descent of the y-curve for that wavelength (figure 3.1) at 
values of & which correspond to raindrop diameters. However, wavelengths 
shorter than 3 cm are not used in rain because of attenuation. At 3 cm wave- 
length, the appropriate values of oO are somewhat lower, and on the basis of 
figure 3.1 we should expect the exact cross-section to be above the Rayleigh 
value for heavy rainfall and below it for light. Haddock's curve in 


figure 4.1 does not do this, but the difference between either of these 


"exact" results and the Rayleigh curve does not exceed 2 db in any case. 
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FIG. 4.1 - Exact and Rayleigh approximate back-scattering cross-sections 
per unit volume of rain-filled space, plotted against rainfall 
intensity. 
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The curves of figure 4.1 are for a temperature of 18C. In the Rayleigh 
region, temperature change will not affect the scattering greatly, since the 
quantity \K{2 to which the Redtearing is proportional does not vary by more 
than a few percent. Out of the Rayleigh region, the scattering is a more 
complex function of the refractive index, but according to Ryde (1946) the 
temperature effect on scattering is still not more than a few percent, 

Several workers have performed experiments to confirm equation Cra 
Marshall, Langille and Palmer (1947) obtained proportionality between P,. 
and Z at 10 cm wavelength, but did not measure the absolute power level. 
Hooper and Kippax (1950) confirmed the wavelength dependence at 10 cm, 3 cm 
and 1.25 cm, and the proportionality to pulse length h. They also reported 
that the absolute value of received power agreed with the theory. Hood (1950) 
found the expected proportionality between P. and Z, but Pe was 3 or 
lL, db too low at 10 cm and 7 db too low at 3 cm. He attributed this to a 
fixed error in the radar calibration or in the method of measuring Z. However, 
Austin and Williams (1951) calibrated a radar set with a known spherical 
target, thus eliminating the need for power measurements, and still obtained 
results 5 to 7 db too low on rain at 10 cm. They also ascertained that 
Hooper and Kippax were measuring peak received powers instead of the average 
power PL, which must have been several decibels lower. Each of these 
experiments show considerable scatter of data within itself, so the experiments 
cannot be said to disagree significantly with each other. There is no doubt 
that they disagree with the theory, and no satisfactory explanation has been 


reached. 


Zove 


4.3 Attenuation by Rain 
For rain, the Rayleigh approximation for attenuation does not apply 


except for wavelengths of 10 cm and longer. Hence the attenuation through 
rain is not given by equation (1.6), but depends on diameter as well as 
water content. 

The data published by Laws and Parsons has been used by Ryde (1945, 
1946) and by Haddock to predict the attenuation at various wavelengths by a 
given rate of rainfall. The authors have repeated some of Haddock's cal- 
culations, using Laws and Parsons! data as plotted by Rigby and Marshall 
(1952) and Haddock's ratio in figure 3.1. In figure 4.2 these various 


calculated attenuations are plotted against wavelength for several rates of 


rai 


FIG. 4.2 - Attenuation versus wavelength 


R for rainfall of various rates, 
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_ Some of the same results are replotted in figure 4.3. Here the abscissa 
is rainfall rate in mm/hr and the ordinate is attenuation divided by rainfall 
rate. The values given by the Rayleigh approximation (equation 1.6) are shown 
as broken lines for comparison. The ordinate in this diagram is the factor k 
in Ryde's (1946) section 7; he suggests that, for many purposes, it may be 
taken as constant for a given wavelength. At 10 cm wavelength, data is scarce 
because Lowan's tables cannot be used, but Ryde's two points, together with 


the knowledge that the Rayleigh value can be used at low rainfall rates, show 
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FIG. 4.3 - Attenuation/rainfall versus rainfall for various wavelengths. 
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that k is very nearly constant, but has such a small value that it is neg- 
ligible for radar purposes. At 3 cm, k is by no means constant and can be 
represented reasonably well by the empirical relation k = 0.009 R03 | where 
R is the rainfall rate in mm/hr. However, if the rainfall rate is brsie to 
lie in a restricted range on any particular occasion, a suitable mean value for 
k could be used. At 9 mm, k is very nearly constant again, though it is 5 to 
10 times the corresponding Rayleigh value. This is because Haddock's ratio is 
near its maximim for these drop sizes at this wavelength. 

Whenever k can be taken as constant, the attenuation in db over the 
path is simply k times the integral of rainfall rate along the path in 
km-mm/hr. 

The data of figure 4.3 is only plotted at widely-spaced wavelengths, so 
it is not suitable as it stands as a design chart for selecting a wavelength for 
a particular purpose. It is found, however, that in figure 4.2, which is a log- 
linear plot, the Ryde curves are nearly straight lines from 3 cm to 10 cm wave- 
length. At 10 mm/hr, for example, attenuation is halved for an increase of 
l.l-cm. At 100 mm/h, the curve is not so linear but the average slope is similar. 

These facts make it possible to interpolate on figure 4.3 and hence predict 
the wavelength at which a given total attenuation will be encountered over a 
given path. In Table 7.1, the data for attenuation by rain for four wavelengths 
at 18C is summarized. 

To determine exactly the effect of temperature on the curve of figures 
4.2 and 4.3, it is necessary to have a set of curves of Q Mie/Q, Rayleigh 
similar to those of figure 3.1 for the refractive indices of water at the 
required wavelengths at each Pobberitanee Ryde apparently computed these 


* This involves lengthy computations, yielding a set of tables the size of 


Lowan's (1949) for each temperature. 


hin 


from first principles for a number of temperatures, since he quotes values 
for the attenuation at other than 18C {Ryde 1945). At that time he used the 
refractive index values quoted by Saxton (1946), some of which differ from 
the more recent values of Saxton and Lane (1952). 

An exact evaluation being too lengthy, the authors obtained approximate 
values for the attenuation at OC for two rainfall rates at wavelengths of 
0.9, 1.24, 1.8 and 3.2 cm, by interpolating amongst the curves of figure 3.1. 
Interpolation is possible since the position of the rising part of these 
curves (the important part for raindrop sizes at these wavelengths) was found 
to depend mostly on the real part of the refractive index. The values for 
attenuation by rain at OC obtained in this way were compared with those at 
18C and plotted as a correction curve in figure 4.4; the temperature correction 
curve of Ryde (1946, figure 7) is shown for comparison. They indicate a 


decrease in attenuation in this wavelength region of about 10% as the tem- 


° McGill estimates 
—— Ryde (1946 - Fig.7) 


Atten. OC 
Atten. 18C 


7IG. 4.4 - Approximate temperature correction curve for attenuation 
by rain (Fig. 7 of Ryde (1946)) with four McGill values 
superimposed. -25- 


perature drops from 18 to O°C, where Ryde found a 30% decrease. However, 
since the McGill values are based on an Miterpodatien whose validity 
could be questioned, it appears that Ryde's (1945) data, from which figure 
4,4 was drawn, must be considered the most reliable available, 
4,4, Shaped Raindrops 
Experiments by Blanchard! on water drops in the laboratory sugyvest 
that in falling the larger raindrops are flattened into roughly oblate 
spheroidal shape with axis vertical. Blanchard also shows that when a 
large drop catches up and collides with a smaller one it will oscillate 
about an equilibrium shape. Such oscillations can also occur as a result 
of aerodynamic forces on the drop as it falls. At any time, therefore, a 
fraction of the total number of drops in the contributing region is non-= 
spherical. Mie's theroy does not apply to non=-spherical particles, and 
no complete theory for spheroids is available. An approximate theory for 
small spheroids which is analogous to the Rayleigh approximations has 
been used by Atlas, Kerker and Hitschfeld (1953) however, and they find 
the the flattening diminishes scattering and attenuation with vertical 
polarization and enhances it with horizontal polarization. There is not 
much data available so far on the shape of falling rain. 
This effect cannot account for the difference between measured and 
predicted average received power because the observations were made with 
horizontally polarized radars, which would lead to greater measured powers 


rather than smaller, 


1 

Blanchard, D.C., 1948: Observations on the behaviour of water drops at 
terminal velocity in air. General Electric Research Laboratories: Project 
Cirrus, Occasional Report No. 7. 


2 
Jones, D.M.A., 1952: Raindrop camera and some preliminary results. 
Proceedings of the Third Radar Weather Conference, McGill University, 
page E-9. 
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5 SNOW 


5.1 Scattering by Snow 

Dry snow consists of ice particles, either single or aggregated. In 
the centimetre region, ice has a dielectric constant of 3.17 which is much 
less than that of water and is independent of wavelength; in the Rayleigh 
region, scattering is proportional to {x and so the small absorption 
coefficient kK of ice can be neglected. It turns out that the back-scatter 
cross-section Y for an ice sphere is only about 0.22 times that for a water 
sphere of the same mass at 3 cm wavelength. 

Although an aggregate flake consists of a mixture of air and ice of low 
mean density, it can be shown (Marshall and Gunn, 1952) that g is the same 
as that for the same mass of ice in compact form. Calculations by Atlas, 
Kerker and Hitschfeld (1953) and by Labrum* indicate that the shape of ice 
particles or snowflakes is unimportant and that the cross-sections are only. 
slightly greater than those for spheres of the same mass (never more than 
twice as great). This means that AT for snow can be predicted by catching 
a sample of snowflakes, melting ane measuring J, po for the resulting drops 
and allowing for the velocity of fall and the per factor, __ 

This was done in Canada by Langille and Thain (1951) who also made radar 
measurements at the same time. The agreement with the form of equation (1.1) 
was within the scatter of the observations, but the received power was 4 db 
less than predicted. As in the rain measurements, this difference is unexplainec 


Langille and Thain's results, on further analysis (Marshall and Gunn, 1952) 


yield an empirical relation between Z and R, where Z refers to unit volume 


- Reported by A. J. Higgs: Proc. Third Radar Weather Conf., McGill University, 
1952, page F-29. 


Soe 


of snowsterm and is the sum of the sixth powers of the diameters of the flakes 
when melted, and R is the snowfall rate at the ground measured in mai aeeae 
of water per hour. The overall relation is Z = 200 Ri-6 identical with 
Canadian results for rain. Since the reflectivity of ice is 0.22 times that 

of water, radar echoes from snow are 0.22 times those from rain of the same 


rainfall rate. 


5.2 Attenuation by Snow 
Attenuation by snow will be partly due to scattering out and partly due 


to absorption of the incident energy (section 3.2). 


According to equations (1.5) and (3.8), 


attenuation (db/km) = 0.4343 2 Qa 
1 
= Oohakas( 5.10.4, 20.) 
spa he "a 


We can use the approximate expressions (3.6 and 3.7) if the wavelength is 
long compared to the particle diameters; the values of (Kj? and Im(-K) to 
be used are listed in table 2.2 and are based on Cumming's (1952) values. 
This gives 


attenuation at OC (db/km) = 0.00349 a + 0.00224 & x Ceca) 
x 


Ryde (1946) suggests that attenuation by snow in the Rayleigh region and 

at centimetric wavelengths is given by the second term only. In Table 7.1, values 
of the attenuation calculated from equation (5.1) for 3 temperatures show that 
the first term (scattering out) becomes important on going to wavelengths 

1 cm and lower. Since Ryde (1946) used refractive index values of Dunsmuir 

and Lamb (1945) his formula for attenuation by snow must be multiplied by 


roughly 2 to get the second term in equation (5.1). 


5.3 Wet Snow and the "Bright Band" 


Atlas, Kerker and Hitschfeld (1953) show that randomly oriented 
small water spheroids scatter and attenuate more strongly than water spheres 
of the same volume; the increase is by a factor of 20 or more if the axes 
of the spheroids are 10 times their diameters. 

Labrum* has computed the scattering by an ice ellipsoid coated with 
water for wavelengths large compared with the size of the particle (corres- 
ponding to the Rayleigh approximation for spheres). He finds that even for 
comparatively thin coatings the composite particle scatters nearly as well as 
an all-water particle of the same shape. He confirmed this result experi- 
mentally. 

Extrapolating from the results of Labrum for small water-coated ice 
spheroids, and Langleben and Gunn®*(1952) for larger water-coated spheres and 
the above-mentioned shape effects for water, it seems not unreasonable to 
assume that the scattering and attenuation by water-coated ice spheroids of 
the size of wet snowflakes at radar wavelengths is similar in magnitude to 
that of spheroidal water drops of the same shape and size. 

This can account for the bright band which always appears in radar pic- 
tures of steady rain. The bright band is due to a region of enhanced re- 
flectivity a few hundred feet thick just below the 0° isotherm. In the light 
of the foregoing discussion, the mechanism would be as follows: 

Snow falling through the freezing level will change from flat or 


needle-shaped particles or aggregates which scatter feebly, to similarly- 


* Loc. cit. 
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shaped particles, which, due to a water coating, scatter relatively strongly. 
As melting proceeds, they lose their extreme shapes and the shape factor 
decreases to 1. Also, falling with greater velocity, the number of particles 
in unit volume becomes less. 

Attenuation, like scattering, is proportional to the number of 
particles per unit volume, and therefore depends on the velocity of fall.. 
Water spheroids absorb more strongly than water spheres, and coated ice 
spheroids presumably absorb like water spheroids. Therefore, absorption by 
melting snow will be many times that of the resulting rain. Although the 
bright band is comparatively thin, serious attenuation my take place when 


radar observations are made through it at shallow elevations. 


6, CLOUDS 


Most cloud droplets do not exceed 100 diameter, so that at centimetre 
and millimetre wavelengths the Rayleigh formulae (3.5) (3.6) and (3.7) apply 


for droplets and equations (1.3) and (1.6) for cloud. 


6,1 Scattering from Clouds 

Cloud droplets being about one-hundredth the diameter of raindrops, 
Ly D® in cloud is about 10° less than in rain. However, equation (1.3) for 
ne received power has nA in the denominator, so that it is possible to 
overcome some of the effect of the very small diameters by going to the 


shortest practicable wavelength. Radar sets working on wavelengths of 1.25 cm 


and around 9 mm are now being used for cloud observations* By directing the 


& 
See "Proceedings of the Third Radar Weather Conference", McGill University, 


Montreal, session B, 
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antenna vertically, the range r, is made a minimum, and most clouds visible 
to the eye are detected by the receiver. 

Radar measures the quantity Z = os po, whereas a more significant 
quantity meteorologically is the Braet content M = z p é D?, To 
relate the two quantities exactly would require the complete drop size dis- 
tribution. However, Atlas and penchant after studying a number of such 
distributions actually observed in clouds, report that a knowledge of the 
median diameter only is sufficient for a fairly accurate relationship, since 
the shape of the distribution is similar for all clouds. Even without this 
knowledge, an approximate correlation between Z and M _ can be obtained 
because the median diameter shows a definite trend with M, and this can be 


incorporated in the empirical law 


Z = 0.0292 uie8*, 


where M is in gm m™3, and Z in mm® m™>, 


6,2 Attenuation by Clouds 
Since the cloud particles lie well within the Rayleigh region at all wave- 


lengths in the centimetre band, equation (1.6) applies. The liquid water 


content M a7 -¥ D, and equation (1.6) becomes 
uf 
attenuation (db/km) = 0.4343 2 M Im (-K) 


for one-way transmission. Thus attenuation by cloud is dependent on the 

liquid water content and is independent of the particle size-distribution. 
Values of the attenuation from water and ice clouds are given in 

Table 7.1 as a function of M for four wavelengths and three temperatures. 


For ice clouds as well as water clouds, Q, << Q,. 


~ loc. cit. pages B 1-8. 
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7. THE ATMOSPHERE 


7.1 Scattering 


Booker and Gordon (1950) have shown that a turbulent atmosphere scatters 
radio waves in a manner depending on the scale of the turbulence, which may 
vary from a few em to a few metres. In particular, long-range point-to-point 
transmission at centimetre wavelengths can take place by such scattering. In 
a typical point-to-point case, however, Goran has shown that rain will con- 
tribute more signal thari the turbulence below about 1 metre wavelength and 
cloud below about 3 cm, provided they fill the space common to the transmitting 
and receiving beams. Attenuation by the same rain or cloud throughout the 
path may however reduce the signal below that for clear turbulent air. 

Booker and Gordon's equation (2v) gives scattered power in any direction. 
It depends on the scale of turbulence A which can be thought of as the 
diameter of the "blobs" of inhomogeneity in the turbulent atmosphere. Making 
the approximation it L/r>>1, substituting the appropriate angle of scattering 
and bringing the definition of ZT into line with section 1 of the present 


1 
paper, we obtain 


2 
2 Fie faye x 107° em@ m3, 
a 


where (AN)2 is the mean square deviation of refractive index in parts per 
million 


and L is the scale in cm. 


k 
Gordon, W. E., 1952: A comparison of radio scattering by precipitation and 
by a turbulent atmosphere. Proceedings of the Third Radar Weather Conference, 
McGill University, page F-17. 
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722 Atmospheric Attenuation 
Van Vleck (1947) has discussed in detail the mechanism of attenuation 


in the clear atmosphere. 

It has been shown by infrared spectroscopy that oxygen molecules 
have a magnetic dipole moment, and water Oy aes molecules have an electric 
dipole moment. These give rise to lines in their absorption spectra 
which lie in the centimetre and millimetre region. This means that radio 
or radar propagation is subject to absorption by the oxygen and water 
vapour in the atmosphere, 

These absorption lines are subject to the usuel emount of "pressure 
broadening". In weter vapour, under atmospheric conditions, the most 
important line has a line half-width of about 0.1 em~l- Since the wave- 
number of the line itself is only 0.74 cn=_, however, the line is relatively 
very broad, having a resonant wavelength of 1.348 cm with half-width of about 
0.18 cm, so that its influence extends over much of the microwave region. 
Additional effects are found due to other absorption lines in the neighbor- 
hood of 0,2 cm wavelength, 

Becker and Autler (1946) measured the damping of a large cavity 
resonator filled with air and water vapour. They used wavelengths between 
0.75 cm and 1.7 cm. Their results are plotted as part of curve (c) in 
figure 7.1. The experiments were performed at 445C, but have been corrected 
to 20C in plotting the curve to bring it into line with the other curves. 
‘The dashed portions of curve (c) are calculations by Van Vleck from a 
theoretical formula with the constants chosen to fit Becker and Autler's 
results. Curve (c) is drawn for one gm m-3, and attenuation is approximately 
proportional to water vapour content in the amounts encountered in the 


atmosphere, 


aioe 


Dicke, Beringer, Kyhl and Vane (1946) used a radiometer to investigate 
absorption by atmospheric water vapour at 1.0 cm, 1.25 cm and 1.5 cm. Though 
the accuracy is not high, the eomniee are consistent with those of Becker and 
Autler. 

The situation for oxygen differs in two respects from that for water 
vapour; the most important absorption is caused by a number of closely~spaced 
lines instead of a single line, and, in addition, oxygen has a line at zero 
frequency, while water vapour has not. 

Van Vleck (1947a) predicted the positions of about 25 oxygen absorption 
lines between 0.45 and 0.56 cm, and these have been measured with a sensitive 
microwave spectrometer at very low pressures by Burkhalter, Anderson, Smith 
and Gordy (1950). At atmospheric pressure, however, pressure broadening converts 
them into an unresolved band of absorption. Lamont (1948) measured attenuation 
in air from 0.45 to 0.58 cm wavelength over a path of a few km in the field, 
and his results are plotted as points in figure 7.1. 

At the time that Van Vleck's calculations were made, it was believed that 
the line breadth constant Ay was the same from line to line, though its 
value was not known. Accordingly, he assigned several alternative values to 
Ay and the predicted attenuation for Ay = 0,02 cml and AY = 0.05 cm 
are plotted against wavelength in figure 7.1, curves (a) and (b). (These are 
similar to Ryde's (1946) figure 1, curves (a) and (b).) 

Van Vleck's predictions include a line at zero frequency, and this accounts 
for most of the attenuation at wavelengths from 1 cm cea ean Its intensity 


is directly proportional to the line breadth and therefore to the pressure, but 


An alternative approach to this kind of absorption can be made by using 
Debye's theory of polar molecules in dielectrics. 


3h 


no measurements appear to have been made to determine the line breadth constant. 
One method is to plot the theoretical curves, choose the one which fits the 
measured values of attenuation in the 0.5 cm region, and assume the remainder 
of the curve is also valid. This makes Av = 0.02 em™/ the best choice 

(curve (a) of figure 7.1) because it fits Lamont's results the best. 

Recent measurements on pure oxygen, however, (see Anderson, Smith and 
Gordy, 1952) show that the line breadth constant varies from one line to 
another in the 0.5 cm region; its value lies between 0.03 and 0.05. We cannot 
say which, if any, of these values applies to the zero-frequency line, so the 
whole attenuation curve for dry air from 1 cm wavelength upward is unknown. 

No measurement has been made above 0.6 cm specifically to define the curve, 
though a single reliable measurement at a wavelength between 1 and 10 cm would 
suffice. | 

The possibility remains of using the results of the water vapour experi- 
ments already mentioned, by extrapolating to zero water vapour content. 
Because of their particular experimental method, Becker and Autler found it 
necessary to eliminate the contribution of oxygen to the absorption in such a 
way that it cannot be determined. The Radiometer method of Dicke et al (1946) 
measured combined attenuation of oxygen and various amounts of water vapour, 
so that if the scatter of the observations were small enough the part due to 
oxygen could be found by extrapolating to zero water vapour. Only the results 
at 1.0 cm wavelength are good enough for this; they indicate a most probable 
value which lies on curve (a) (corresponding to AV = 0,02 cm!) and that 
the attenuation is unlikely to be as great as shown by curve (b) which can 


be taken as an upper limit. 
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FIG. 7.1 - Attenuation by air and water vapour at 200 and 
atmospheric pressure. 
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Temperature and pressure are discussed by Van Vleck in his two papers. 
For a given quantity of water vapour, attenuation is proportional to pul 
and to T-2 10-278/T at the maximum of the resonance curve, to p and to 
7-3 107278/T on the sides of the curves, and to p and to ¢73/2 well away 
from the resonance. Of course, if the air is kept saturated with water vapour 
while the temperature changes, the water vapour content itself will alter, 
and this additional factor must be taken into account. 

Pressure effect on atmospheric oxygen arises from increased pressure 
broadening and also from an increase in the number of molecules per unit 
volume. In the region of strong absorption, the attenuation is proportional 
to p and to 7-3/2, but in the region of the "zero-frequency" absorption 
(from say 2 to 10 cm wavelength) it is proportional to p“ and to q75/2, 

Table 7.1 is a summary of the results of the preceding sections and gives 


working values of the attenuation by rain, snow and clouds. 


on3'Jue 


SNOW* 


R in 

mm hr 
of melted 
water 


TABLE 7el 


363x107 RE 06 
+ 68,6x10" R 


3.3x107 Rte® 
+ 22,9x107°R 


3.3xL07R | 
+ 15.7x10 ~ R 


0.0483 M 
0.0630 M 
0.0858 M 


0.112 M 
(extrapolated) 


2,46x1072 M 
8219x1074 M 


5 963x10"4 


3632x107 4R1 6 
+ 12,2x1074 R 


3432x1074R1-6 


+ 4e06x107% R | 


3,32x1074pL 06 


+ 2,80x1074 R 


436x107? M 
1246x102 M 


1x1073 M 


Working Values of ATTENUATION in db.km~ (one way) 


1.48x1073n1 +6 
+ 1,78x10"? R 


148x107 3pi6 
+ 0.59x1072 R 


1,.48x1073R2 °° 
+ Os4lx1072 R 


6635x107? M 
2.11x1072 M 
1245x1072 M 


5933x107 5RL 26 
+ 244x107 R 


5 633x107 3R2 «6 
+ 0,81x1073 R 


53x07 pled 
+ 0.56x1072 R 


8, 74x107> M 
2.91x1073 M 
2x107? M 


& These are empirical equations. For the more exact relations on which they are based, 
see section 4.3 and figure 4.3. 

kkThe effect of temperature is discussed in section 4.3 and summarized in figure 4,4. The 
effect of non-sphericity on Q. and Q, has been neglected. b 

* As long as snowflakes are in fhe Rayleigh region. A value of R#10 or Rie 
limit for snowfall rates. 


=39 is an upper 
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8.1 


8.2 


8. THE DIRECTION OF FUTURE RESEARCH 


Refractive Indices 


The measured values of the refractive index of water for various 
wavelengths and temperatures are quite adequate for microwave- 
precipitation work. The fact that the measured values fit the 
theoretical formulae so well (see, for example Saxton, 1946) allows 
interpolation and extrapolation about these values to be made with 


confidence. 


The situation for ice is not so well established. There is 
agreement on the value of the real part of the refractive index and 
on its independence of temperature and wavelength variation in the 
microwave region. However, measured values of the imaginary part 
differ by a factor 2. Probably further determinations of this 


quantity are in order. 


Rain 

While anomalies in raindrop 8ize-distributions continue to be 
reported, they seem to appear at an equilibrium rate and make a very 
small dilution of the general broadagreement. Data will cmtinue to 
accumulate as an adjunct of any scattering and attenuation measure- 
ments that are made. For a better understanding of the rain mechanism, 
a more detailed study of the size-distribution with time from a given 


shower, and variations with precipitation type, is needed. 


Though Ryde (1945) evaluated the attenuation by rain at four 
wavelengths, and various temperatures, the work may bear repeating for 
the slightly changed values of the refractive index of water which 


have been reported since his computations were made. The final results 


-39= 


8.3 


should not be very different but in view of the appreciable 


attenuation by rain, and the magnitude of the computational work, 


ait might be worth while essentially repeating Ryde's work by 


compiling a table of the Mie amplitude functions for water at OC 


similar to Lowan's (1949) 18C table. 


If non-spherical drops occur in natural rain, the scattering 
and attenuation could be considerably affected. Shaped raindrops 
are thought to be rare but a raindrop photography program undertaken 
recently by the Illinois State Water Survey should provide a definite 


answer to this problem. 


There is a general scarcity of point to point measurements of 


attenuation by rain. 


One outstanding unsolved problem is the source of the 4=7 db. 
difference between the actual signal received from both rain and 


snow and that predicted from the theory. 


Snow 

The dielectric properties of ice and snow and wet snow have 
been measured in the lab by Cumming (1952). There is still, however, 
the disagreement about the value of the imaginary part (k) of the 
refractive index. Labrum in Australia has observed the scattering 


from single melting spheroidal ice particles. 


There has been little experimental work done on either scattering 
or attenuation by falling snow. There is a need for laboratory 
experiments to study the melting of natural snowflakes, both single 
crystals and aggregates. Point to point measurements of attenuation — 


by melting snow would be particularly valuable. 


=)0- 


8.4 Cloud 


8.5 


Laboratory experiments on the scattering by single water 
spheres have been done, and give results in good agreement with 
theory. Laboratory experiments to determine the dielectric 
properties of water cloud are possible and might provide amongst 
other things a possible clue to the 4-7 db discrepancy between 


actual and predicted scattering mentioned above. 


Atmosphere 


The attenuation due to water vapour has only been measured 
at wavelengths 0.7 and 1.7 cm, but the necessary theoretical basis 


for extrapolating beyond these limits seems to be satisfactory. 


The most noticeable gap in present knowledge of atmospheric 
attenuation is a direct and reliable measurement of the zero-frequency 


line in oxygen under atmospheric conditions. 


A laboratory experiment at any wavelength over 1 cm on dry air 
at atmospheric and at very low pressures, or a field experiment at 
any wavelength over 3 cm at a low humidity would determine the 
line breadth constant and remove the uncertainty at all wavelengths 
above the 0.5 cm band. However, since the attenuation appears to be 
not much greater than 0.01 db/km it is difficult to measure and at the 
Same time usually unimportant in practice. Any application in which 
it becomes important to know it precisely will probably suggest a way 


of measuring it accurately. 
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INTRODUCTION 

Part I of this report has, among other things, provided the basic 
pormclae from which attenuation coefficients may be computed, tn addition, 
it has included many of the most recent computations of these coefficients, 
It is now our purpose to determine the cumulative effects of the various 
atmospheric constituents on microwave attenuation, Here our attention will 
be restricted to the very short wavelengths (“. 3.2cm). 

As indicated in the preface we shall compute the losses to be expected 
under a number of model atmospheric conditions ranging from a clear standard 
atmosphere to a severe line squall. xcept in the case involving the line 
squall, the models are simplified in that horizontal stratification and 
homogeneity of the various cloud and precipitation layers are assumed over 
the area covered by the radar. The assumption of horizontal homogeneity, 
particularly of precipitation, is generally invalid in reality. However, 
-the complexity of the precipitation distribution is generally such as to 
make this simplification necessary in order to arrive at a satisfactory 
estimate of the attenuation, The resulting values should therefore be 
considered as typical of that to be expected in a particular case. 

The model weather situations which we have chosen are based primarily 
on physical considerations and are therefore largely independent of gzeo= 
graphical location. The frequency of occurrence of each type of situation 
varies widely from place to place, and we have made no attempt to treat such 
climatological problems, Where we have referred to cloud densities or rates 
of precipitation as "high, low, or moderate" the adjective applies to con= 
ditions in temperate United States, unless specifically stated otherwise, 

Although we have selected what are believed to be reasonable models, the 
validity of.some of their physical details may be subject to question, For 


example, sharp boundaries have been assumed between all layers, and the cloud 


density or precipitation rate is regarded as constant throughout each layer. 
These simplifications are required for computational facility and have little 
effect on the final results. A somewhat more. critical assumption is the 
existence of substantial amounts of supercooled water in the presence of ice 
crystal cloud, There is some controversy regarding the validity of this 
supposition, although there is little argument that such conditions must exist 
at least in the early stages of significant. temperate latitude nratigin scene 
In any case, the effect of modifying such premises may be estimated from the 
data provided. 

2. BISCUSSION AND RESULTS 
2.1 _ ase 1 


._ Standard (Maritime Polar) Air; Vase IT, Maritime Jropical Air 


In order to Sasi sk a base-of reference, we shall start with the U.5, 
Standard Atmosphere (Beers, 1945). Inasmuch as this refers to a dry atmosphere, 
we have added water vapor at each level in accordance with the average of a 
number of typical air mass soundings. These soundings included all the con— 
ventional air masses other than continental tropical (cT) for both winter 
and summer, The resulting vertical temperature and water vapor distribution 
which is tabulated at the left of Figure 2,3 resembles ad aubeen eee ie 
polar (mP) air mass. This is Case I, 

Case II. covers a maritime tropical (mT) air mass whose temperature and 
water content is tabulated ati the left of Figure 2,4, ‘the temperature data 
was obtained. from Figure 18 of "The Thunderstorm® (1949) and represents the 
typical temperatures to be expected surrounding a middle latitude thunder= 
storm, Since no water content values were specified in the figure, a water 
content lapse rate for typical summer maritime Gulf air was used. This water 
vapor distribution was obtained from a study by Plank (1952, Figure 13) and 


is based on a consideration of some 200 soundings. 


In the attenuation computations for both Cases I and IT as well as in 
all those cases which follow, we have attempted to use the best available 
estimates of the attenuation coefficients for oxygen and water vapor. However, 
in ee the literature, one finds considerable discrepancy between the 
various determinations of these coefficients. 

This problem has been discussed to some extent in Section 7.2 of Part 1,* 
It was noted there, for example, that two estimates are available for the at= 
tenuation due to oxygen. Although the lower of these appears to be slightly 
favored by recent evidence, this is not certain. In order to obtain an 
estimate of what is probably an upper limit for the oxygen losses without 
entirely disregarding the evidence favoring the lower coefficient, we have 
split the difference between the two oxygen curves of Figure 7.1 of Part I, 
This is sho-n in Figure 2,1 where the coefficients are plotted in decibels 
per mile for a two-way path, All of our initial computations and graphs 
are based on this intermediate estimate, Later (Section 3.3) we shall 
consider how the results are affected if the lower estimate is used, 

The attenuation coefficients for water vapor have been taken from the 
theoretical values of Yan Vleck (1951, Table 8.2 and Figure 8.6) and are 
shown as a function of wavelength in Figure 2.1, These values are in es= 
sential agreement with the water vapor attenuation values one obtains from 
the Becker and Autler curves in Figure 7.1, Part I, Only at 3.2 cm do 
the values differ appreciably; however, since the water vapor losses at | 
this wavelength are very small, no significant modification of our results 


appears necessary. 


*Also see Appendix, 


Fig. 2e1 


The variation in oxygen and water vapor attenuation with 
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wavelength. 


Fig. 2.2 


The variation in attenuation with altitude in a standard (mP) 
atmosphere and a maritime tropical atmosphere. 


The oxygen curves apply to both atmospheres within about t 10%. 


_ The right hand water vapor curves apply to an mT atmosohere; 
those on the left apply to the standard (mP) atmosphere, 


The 3.2 cm water vapor curve for a winter continental polar 
Similar curves apply to 
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In Figure 202 the two-way losses for Cases I and II have been plotted 
for one-mile increments of altitude assuming the oxygen and vapor densities 
at the middle of the layer to be the mean for that layer. The variation of 
oxygen absorption with altitude is based on a statement by Van Vleck 
(1951, p 647) to the effect that the oxygen coefficient varies as the quo- 
tient of density by temperature.* Inasmuch as the density—temperature 
quotient varies but slightly at any constant level in even radically 
different air masses, the same oxygen curves may be used with satisfactory 
accuracy (within about + 10%) for practically all air masses, 

The differences between the various air masses are due to their 
contrasting vapor content and are-illustrated in Figure 2.2. For comparison, 
dry continental polar (cP) air is also shown there, At the surface, for 
example, the Hj0 vapor coefficient for mT air is mare than-20 times that 
for cP air, However, the difference between the H50 curves for the standard 
(mP) atmosphere and that for the mT air is everywhere less than a factor 
of 1.7. Since the water vapor loss to be expected in cP air is so slight, 
it would be quite appropriate to consider oxygen attenuation alone when 
dealing with such an air mass, This is especially true at 3.2 cm. (The 
continental polar curve in Figure 2.2 is drawn only at a wavelength of 3.2 
cm to avoid confusing the diagram, The effects are proportional at the 
other wavelengths, ) 

Water vapor attemation falls off much more rapdily with altitude than 
does that for oxygen, deereasing by more than one order of magnitude (in 
Cases I and II) from the ground to an altitude of 4 miles.._Also, relative 
to oxygen, water vapor is of about equal importance in the lower 4 miles of 


the atmosphere at wavelengths of 0.9 and 1.8 cm, up to 7 times more important 


*This is evidently incorrect. See Appendix 
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Cumulative attenuation due to a maritime tropical atmosphere. 


at 1.25 cm, and more than 5 times less important at 3.2 cm, depending on 
the air mass... Above 4 miles, oxygen soon takes over, with the water vapor 
effect becoming: small or negligible above about 6 miles at all wavelengths. 

The total two-way attermation accumulated from the ground up to any 
particular altitude is given in Figure 2.3 for Case I, the standard (mP) 
atmosphere, and in Figure 2.4 for Case II, the mI atmosphere, The totals 
do not exceed 0.8 decibels except at 1.25 cm where they are about 1.4 = 
1.8 db. The altitudes at which the cumulative losses (upwards). are equal 
to 67% and 90% of the 12-mile total are indicated on Figures 2.3 and 2.4. 
Except at 1.25 cm, 67% of the total attemation occurs at an altitude of 
about 3-4 miles and 90% occurs at about 7-8 miles in both cases, Due to 
the greater. importance of water vapor at 1.25 cm, the 67% and 90% altitudes 
are lowered considerably. It is to be noted that these air masses are 
more opaque to 1.25 cm than to 0.9 cm. The same is true for any cloudless 
air mass having a significant amount of water vapor. 

The differences in total attenuation that would result. in Cases I and 
It if Van Wleck's lower estimate for oxygen were used: to eeclacats inter— 
mediate estimates are shown in Table 3.3, The lower coefficients decrease 
the total losses by 0.053, 0.067, 0.175, and 0,22 db at the four wavelengths 
of interest, the values increasing toward the shorter wavelengths, Table 
3.3 shows that these reductions correspond roughly to 30% of the total loss 
at the two extreme wavelengths, to 17% at 1.8 cm, and eens 10% at 1.25 
em for both Cases. I and II, 
2.2 Non-vertical. Paths 

In looking vertically downward, total attenuation is a function of 
altitude. Along:a path other ae vertical, the total loss from a given 


altitude in a.horizontally stratified atmosphere is magnified according to 
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Correction factors to be applied to any of the cumulative 


attenuation figures to convert their values from vertical 


incidence to values applicable to observation at horizontal range 


R and altitude h, 


100 


the ratio of the oblique to the vertical path lengths or by the secant of 
the intervening angles. On this basis, it is possible to draw a series of 
eraphs giving a factor C by which the vertical attenuation from any given 
altitude is increased if the radar is directed at a ground point a horizon- 
tal distance of R miles away. Such a series of correction curves (computed 
for a flat earth) is shown in Figure 2.5 for aircraft altitudes ranging 
from 2 to 12 miles, It might be noted here that these curves are equally 
applicable to ground radars directed upwards toward airborne targets. To 
illustrate the use of this figure, Case I shows the total vertical two-vway 
attenuation from § miles to the surface at 0.9 em to be 0.66 db. If the 
radar were directed at a ground point 100 miles (horizontally) distant the 
total loss would be increased by a factor of 12.6 to a value of 8.3 db. 


Errors in Figure 2.5 due to the flat earth approximation are less than 0.1%. 


2.3 Case III- An Ice Cloud_in Maritime Polar Air Overrunning Modified 
Continental Polar Air 


In this case, we consider an ice cloud imbedded in standard (mP) air 
above a mixture of continental and maritime polar air. The vertical distribu- 
tion of temperature and absolute humidity is given at the left of Figure 2.6. 
The ice cloud has a density M of 0.2 gm/m? and extends from an altitude of 2 
to 5 miles; conventionally it would be called altrostratus although cirro- 
stratus may well have similar characteristics. An ice content of 0.2 em me? 
is thought to be on the moderate to high side for such clouds. Saturation 
with respect to water (100% relative humidity) is assumed in the cloud layer, 
thereby providing favorable conditions for the growth of the ice particles, 
However, unless new air was being brought into the layer, the supply of 
moisture would soon be depleted and further growth would cease, An ‘ice 
cloud in earlier or more active growth would actually exist in conjunction 
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Cumulative attenuation due to the conditions described in Fig. 2.6. 
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with a supercooled water cloud. 1m that event, the attenuation in the cloud 
deck would be substantially greater. The present case may be considered as 
about the upper limit for an old ice cloud, 

The relative magnitudes and variation of the individual components. of 
attenuation with altitude is shown in Figure 2.6. Oxygen attenuation is 
assumed to be the same as in the previous cases although all the curves should 
actually be. shifted slightly to the right due to the lower temperatures. 
Attemation by water vapor is substantially lower than in Cases I and Il 
because of. the colder and drier air below the cloud. The inversion in the 
H90 vapor attenuation coefficient reflects the increased. moisture content 
in the vicinity of the cloud deck. (No attempt was made to draw the points 
of discontinuity with a high degree of accuracy; the curves are drawn for 
values of vapor density computed at the midpoints of layers one mile in 
thickness.) 

Attemation by the ice crystal cloud was computed on a slightly different 
basis than that indicated in Table 7.1 of Part I, In the latter table, the 
coefficients for ice cloud assume that the particles are all small relative: 
to the wavelength and that there is no loss due to scattering. This is quite 
satisfactory at a wavelength of 3.2 cm and the coefficient shown in Figure 2,6 
for this wavelength agrees with that of Table 7.1 for a temperature of 0°C, 
There is, however, evidence that some ice clouds may be composed of relatively 
large particles, of the order of a few millimeters, These would produce some 
loss due to scattering and should more truly be referred to as "snow clouds", 
In order to obtain an upper estimate of this effect, the ice content was 
converted to an equivalent precipitation rate by means of the relation 
R = 3.6 Mv where R is the rate in mm/hr, M is the ice content density in 


em/m?, and v is the particle fall velocity in m/sec, A value of 1 m/sec was 


LL 


chosen for v.. The attenuation was then computed using the relations listed 
in Table 7,1 (Part I) for snow. The greatest differences occur at the 
shortest wavelength where our method gives an attenuation approximately twice 
that for a small particle ice cloud, Although this approach may overestimate 
the effect of scattering, the increase in attenuation tends to compensate for 
the omission of the shape effects from Table 7.1. 

Another effect which was neglected in our computation of the ice cloud 
losses was..the variation of the attenuation coefficient with temperature. 
If the loss js due primarily to absorption, as it is for small particles, 
then there is indeed a variation with temperature corresponding to the 
ig ra of the complex portion of the refractive index, However, losses 
attributed to scattering by ice are not temperature sensitive since they 
are a function of the real and constant part of the refractive index. 
Throughout this paper we have used the ice attenuation coefficient cor- 
responding to a temperature of 0°C, Neglect of the temperature variation 
again results in an upper limit for attenuation by ice, 

The detailed consideration of the above effects may be of somewhat 
academic interest when it is noted that the contribution by the ice cloud 
to the total loss over the 12-mile vertical path is less than 6% for a : 
wavelength of 0.9 cm and of the order of 2% at the other wavelengths, It 
_is to be emphasized that these figures are the result of high estimates of 
the ice content for such clouds as well as the largest probable attenuation 
coefficients, Indeed, it is exactly for this reason that we may neglect 
the effects of shape and temperature, For most practical purposes, there— 
fore, the attenuations predicted by the formulas of Table 7.1 (Part I) will 


be of sufficient accuracy, 


ire: 


Figure 2.7 shows the variation of the cumulative. attenuation from the 
ground up at the four wavelengths of interest. The absence of a noticeable 
change in slope (except at 1.25 cm) at the height of the ice cloud graphically 
indicates the small losses in this cloud, The noticeable change at 1.25 cm 
is associated with the inversion of vapor content in the cloud and not with 
the ice itself. It is interesting to observe that the magnitudes of the 
cumulative losses are everywhere smaller than in Cases I and II because of 
the lower vapor eerrenes However, the 67% and 90% altitude is only slightly 


higher than in the previous cases. 


2.4 Case IV — Heavy Snow Falling trom a Dense Mixed Water-lce Cloud 
in. an Atmosphere Otherwise Similar to Case III 


The situation in this case is a more active form of Vase III in which 
the "snow" cloud has doubled its density and is growing at the expense 
of a coexisting supercooled water tloud of density 0.1 em/tn? « Snow is falling 
from this mixed cloud at a rate equivalent to 5mm of water per hour, which 
corresponds to a rate of accumulation on the ground of roughly 2.5 inches of 
snow per hour (depending on its density), This is considered to be a heavy 
snowstorm, ‘The temperature and vapor distribution is identical to that of 
Case III, with a ground temperature of O0°C, 

There may be some question as to the reality of the combination of 
physical conditions in the assumed model, In particular, one may ask whether 
such a heavy snowfall may occur under these conditions, he answer depends 
on rather complex computations of the growth of ice crystals by diffusion and 
coalescence with cloud, and upon some nebulous corrections for aggregation 
of the particles with each other, for this type of study such a treatment 
hardly appeared worthwhile, particularly in view of the variability of the 


atmosphere. ‘therefore, this model and those which follow are based on what 
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Cumulative attenuation due to the conditions described in Fig. 2.8. 


we believe to be reasonable estimates according to the few available observa— 
tions (Cunningham and Miller, 1948; Cunningham, 1952) and approximate 
theoretical considerations, 

The relative magnitudes and variaticn of the individual components of 
attemation are shown in Figure 2.8. Oxygen and water vapor losses are 
identical to those of Case LII, Attemation by the ice crystal or snow 
Cloud was computed as previously and is about twice that of Vase III, The 
loss due to the heavy snow below the cloud was computed according to 
Table 7.1 (Part I) for a temperature of 0°C. This loss corresponds to 
that of a snow cloud having a density of 1.4 gm/m? More realistically, 
the ice crystal attenuation should increase downward and join continuously 
with that due to the snow thereby reflecting the growth of the particles 
and increas* in precipitation rate downward. However, the results would 
be substantially unchanged. 

The losses due to the supercooled water cloud were computed by extra- 
polating the temperature variations shown in Table 7.1 (Part I) for wave- 
lengths of 3.2 and 1.24 cm.* The temperature variations at 0.9 and 1.8 cm 
were assumed to be proportional to those at 1.24 cm. Note that the loss in 
liquid cloud increases with decreasing temperature; the opposite is true for 
ice clouds, It is evident from the summary in lable 3.1 that in this case the 
liquid water is by far the major contributor to the attenuation at all wave- 
lengths; it accounts for over 46% of the 12—mile total at 3.2 cm and as much as 
63% at 0.9 cm. 

The losses accumulated upwards are shown in rigure 2.9. (Note the change 
of scale on the abscissa.) ‘ihe rate of change of the cumulative curves is 


greatest in those layers occupied by the most lossy media; i.e. liquid water. 


* The cloud loss values for 1.24 cm were used instead of those for 1.25 cm, 


iMate 


This is most evident at the wavelength most affected by liquid water, namely 
O.9 cm, Because the liquid water loss is greatest at the colder temperatures, 
the slope of the curve is greatest near the cloud top. For this reason, the 
67% and 90% cumulative losses are ccncentrated in that vicinity. 

it should be noted that the absence of the liquid water cloud would 
reduce the cumulative losses by more than 46% at all wavelengths, 1n that 
event, the resulting attemations would be as shown in Table 2.1 where’ thew 
are compared to the Case III totals, ‘the differences, which represent es= 
sentially the losses in the heavy snow, become significant only at the two 
shorter wavelengths, 1[t is evident therefore that losses due to all but 
the heaviest snowfalls are small except at the shortest wavelengths. 
Vertical Path (db> 


Table 2.1=lotal Two-way Attenuation Over 12 mi 


Attenuation 
1 a ES et igus CEE 1.8 1,25 0.9 


Case LV without 
water cloud 0201 0328 oJae 1. OS9 
Gase III eked eh Se Pad 2036 
Difference soy gOS: olhy5 ob bh 
Note also that in this case C.9 em is subject to more loss than 1.25 cm 


at all altitudes in contrast to the previous cases, This is due to the 
increased loss (at 0,9 cm) in the snow in the lower levels, and in the mixed 
clouds, The only advantage which 0.9 cm has over 1.25 om is the smaller 
water vapor losses, In general, it appears that the atmosphere becomes more 
cpaque to 0.9 cm than to 1,25 cm when water vapor contributes less than 
roughly 50% to the total 1.25 cm loss (See Table 3.1.) It will be seen 

that this is the case whenever there is any significant water cloud or rain 


present or snow of greater than moderate intensity. 
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2.5 Case V —_A Thick Dense Water Cloud in Maritime Polar Air 


The model atmosphere assumed in Case V corresponds to the standard 
maritime polar air mass of Case I into which we have inserted a supercooled 
nimbostratus cloud having a liquid water content of 0.2 g/m? throughout its 
3-mile depth, (Saturated conditions within the cloud also require an increase 
in the vapor density over that in the clear standard air.) The term "nimbo- 
stratus" is Madstinies restricted to precipitating cloud; here it is used to 
denote only a thick stratiform deck. An average of 0.2 gm/m throughout the 
cloud is somewhat on the high side although it would be considered about 
average for a somewhat thinner supercooled stratus or stratocumulus cloud. 

The magnitude and vertical variation of the loss through the various 
constituents is shown in Figure 2.10, Here again the liquid cloud contributes 
most to the total loss at 0.9 and1.8 cm, This may more readily be seen in 
the summary of Table 3.1. At 3.2 cm, the loss in cloud is about equal to 
that of oxygen, while it equals the effect of water vapor at 1.25 cm. 

The accumulated losses upwards are shown in Figure 2.11. The attema- 
tion in the cloud alone is roughly equal to that through the entire 12-mile 
thickness of clear standard atmosphere except at 0.9 cm where it is nearly 
three times as great (see Table 3.1). Here again the effect of the liquid 
cloud crowds the 67% and 90% contours together at a level near the cloud top 
at the shorter wavelengths, The decreasing importance of liquid water with 
increasing wavelength is indicated by the greater spread of these contours 
at 3.2 cm. 

Case V is of particular interest in that it represents a situation in 
which the one advantage of 0.9 cm over 1.25 cm (i.e, lower attenuation 


through water vapor) is almost exactly compensated by the greater loss in 
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Cumulative attenuation due to the conditions described in Fig. 2.10. 
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water cloud at the lower wavelength. For clouds of any greater density and 
for rainfall of any appreciable intensity (or snowfall of greater than mod= 
erate intensity), 1.25 cm is superior to 0.9 cm. This is especially true 

in situations where water vapor is of relatively minor importance, as in Case 
IV. Incidentally, it may be noted here that 1.25 cm radar is also less 
severly affected than 0.9 cm radar by scattering from clouds and precipita- 
tion, thereby being freer from weather clutter, 


2.6 Case VI — Heavy Snow Turning 
Cloud, 


In this case, hy pean of Case V has turned partly to snow which is 
falling at a rate of 4 mm/hr and becoming rain at the melting level with a 
rate of 4.5 mm/hr. The snow has dried out the cloud above the 0°C level, thus 
leaving only 0.05 gm/m? of liquid water, Below this level, the cloud still 
retains its density of 0.2 gin/m? « Without replenishment, this cloud would 
soon be "washed out" by coalescence with rain. The situation is probably 
typical, in.a simplified way, of continuous rain of moderate to heavy intensity. 
For example, a rain of 4,5 mm/hr is exceeded about 10 hours per year in Eng- 
land (Bilham, 1938). 

The losses due to each of the elements other than rain were computed as 
described previously. Their variation with height is shown in Figure 2,12, 
(The abscissa vata here have been displaced one cycle to the left in order 
to plot the high rainfall coefficients.) The rainfall coefficients used in this 
section were obtained from Figure 4.2 of Part I for a temperature of 18°C, It 
should be noted that the coefficients obtained from the latter figure may not 
agree exactly with the coefficients one would obtain by using the empirical 
equations given in Table 7.1. (Part I), although these equations are accurate 


to within 10%. 
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Cumulative attenuation due to the conditions described in Fig. 2.12. 
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Also, no corrections were made for variations with temperature, Although 
such corrections have been shown to be relatively small at the wavelengths 
with which we are concerned (See Section 4.3, Part I) the magnitude of the 
corrections alone may exceed the magnitude of the total contributions of the 
less effective atmospheric constituents, Thus, if serious consideration is 
to be given to the quantitative values of the lesses due to these minor elements 
(such as are listed in Table 3.1), the temperature correction for rainfall 
should be included, In considering total losses, however, such corrections 
could not decrease the results by more than about 8% at any wavelength, 

It is evident from Figure 2.12 and from Table 3.1 that rain is by far 
the largest contributor to the total loss, Table 3.1 shows that the rain is 
responsible for more than 50% of the total attenuation at all wavelengths. 
This is indirectly.evident also from the cumulative curves of Figure 2,13 
where 67% of the total 12-mile loss occurs below an altitude of 2 miles, 

It is of interest to note, however, that at 3.2 cm oxygen still causes a 
substantial part of the loss. in other words, it takes some fairly heavy 
water clouds (Cases IV and V) or some moderate rainfall to equal the effect 
of oxygen at 3.2 cm, The cumulative losses become quite substantial in this 
case, particularly at the shorter wavelengths. Even at 3.2 cm, however, a 
loss of almost 0.6 db becomes significant when we consider that itis 
magnified by about 12.6 for an nienrate 8 miles high and a ground target at 
a range of 100 miles, If the ground target were barely detectable ina 
hypothetically transparent atmosphere, such a loss (7.5 db) would require 
an increase in overall system sensitivity of about 5.6 in order for the tar— 
get to remain detectable, Similarly, if one wished to increase the range 


of a 3.2 cm radar from 50 miles to 100 miles (from an altitude of 8 miles) 
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Cumulative attenuation due to the conditions described in Fig. 2.14. 
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under these weather conditions, a fourfold increase in overall system 
sensitivity would be required to compensate for the inverse square range 
effect and an additional factor of about 2.4 to compensate for the increased 
atmospheric losses, or a total increase by almost a factor of 10, The 


corresponding power. increases at the other wavelengths would be even greater. 


2.7 Case VII — Viewing the Ground. from above an Intense Thunderstorm in 
Maritime Tropical Air 


This case depicts one of the most extreme conditions likely to be en= 
countered, it corresponds to a heavy shower or thunderstorm imbedded in a 
moist maritime tropical air mass. Heavy rain is falling at a rate of 25 mm/hr 
from an altitude.of 3 miles, A dense water clond with a content of 0.5 gn/m? 
is present between 1 and 6 miles and a thick ice crystal cloud of 0.4 g/m 
content exists in the tep 2 miles corresponding to the anvil of a cumulonimbus, 
In the computations, liquid water (0.1 gm/m?) was assumed to coexist with the 
ice cloud although with the low temperatures: existing at this altitude such 
an assumption is questionable, In any case, the liquid water near the cloud 
top accounts for less than 5% of the total 12—mile loss at all wavelengths, 
Saturation with respect to water is assumed from the ground to the cloud top, 
The temperatures were taken directly from Figure 18 of "The Thunderstorm" and 
represent typical conditions surrounding a mid=latitude thunderstorm, Absolute 
humidities within the cloud are prescribed by the temperatures and assumed 
saturated conditions, 

Such intense storms very rarely have dimensions greater than 5-10 miles, 
although a series of similar cells may extend for as much as 200 miles or more 
in one direction in the form of a line squall. Thus, unless the radar happens 
to be directed along the larger dimension of such a line, it is not possible 
to consider these storms to be horizontally homogeneous, The case of a radar 
looking at. the. ground. in a direction perpendicular to a line squall is covered 
in Case VIII. 
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Variation of total attenuation with horizontal range for a 
radar viewing the ground through an intense thunderstorm in a 
manner shown by the sketch at the top. 


The results for Case VII are shown in Figures 2.14 and 2.15. The relative 
importance of the panidas constituents in contributing to the 12-mile total 
loss is shown in Table 3.1. It is quite evident that such intense rain and dense 
Liquid cloud overwhelm the effects of oxygen, water vapor, and ice or snow, 
It should be noted that these attemations apply only along the vertical 
path and should not be increased by the correction factors of Figure 2.5 
unless one is interested in the magnitude of the losses looking through 
the long dimension of a line squall. In this direction, the total two- 
way loss from a height of 8 miles out to a ground range of 100 miles 


would be about 80 db even for 3.2 cm, 


2.8 Case VIII -. Views 

A more likely situation than that of Case VII would be one in which 
the radar is viewing the ground from a position behind a thunderstorm, 
For this case, we assume the same vertical distribution of the various 
constituents as in the previous one. However, the line squall is given 
a horizontal thickness of 10 miles in the direction of the radar beam, 
The aircraft is at an altitude of 12 miles and is located 10 miles behind 
the storm. 

The total loss out to any point on the ground depends on the part— 
icular path through the storm as indicated by the sketch at the top of 
Figure 2.16, Computations were carried out for ground ranges up to 50 
miles from the position of the radar and the results are shown in the 
lower portion of Figure 2.16. For ground ranges less than 10 miles the 
beam is entirely in the clear maritime tropical air, Beyond 10 miles, 
the path through rain increases rapidly out to 13 miles with a corre— 
sponding rapid increase in the loss; between 13 and 20 miles the distance 


in rain increases slowly while that in cloud is increasing rapidly thus 
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decreasing the slope of the loss curves; at 20 miles the distance through 
rain is a.maximum and that through cloud nearly a maximum resulting in 
the maximum losses, Similar discontinuities in slope are evident beyond 
the peak, The maximum losses are considerable, If the target was under 
or just beyond the storm, it would be blocked quite effectively, possibly 
even at 3.2 cm, However, since the storms would themselves be displayed 
on the radar scope, there is the possibility of vectoring the aircraft 
through them to a more pence position, 

It mast be emphasized, however, that an airborne radar operating in 
the vicinity of a complex distribution of intense cloud or precipitation 
cells would be subject to heavy and variable attenuation as the beam 
scanned, A ground radar pattern viewed under such conditions would con- 
tain shadews of variable intensity behind each of the storm cells, in 
addition to the echo clutter of the storms themselves, The resulting 
display could be extremely confusing, if not entirely worthless, 

2.9 Approximate. Statistical Data | 

Heretofore we have intentionally neglected to consider the freq= 
uencies with which the various model weather situations could be expected 
to eccur, Such considerations are necessarily dependent on climatologi- 
cal studies pertinent to a particular geographical region. Furthermore, 
these studies must incorporate data on the appropriate atmospheric para=- 
meters throughout the depth of the navigable atmosphere as well as at 
the ground, Kwen the climatological information on surface rainfall 
intensities is inadequate, not to mention the absence of satisfactory 
statistical data. on cloud density and thickness, 

Nevertheless, some rough estimate of attenuation expectancy must be 


made in order to furnish some basis for evaluation of the severity of 
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the problem. For this purpose we have utilized Bilham's (1938) data 

on surface rainfall intensities and duration in England, with the full 
realization that attenuation statistics based on this data can be repre- 
sentative of but a restricted portion of the globe. These, particular 
data have been used extensively in studies similar to the present one, 
notable by Ryde* (1946) and more recently by Rogers* (1953). 

In estimated attenuation through the vertical extent of the 
atmosphere from surface rainfall data, we shall make two very rough.ap- 
proximations: (1) that the general characteristics of the precipitation 
model of Case VI (Figure 2.12) are typical of a wide range of rainfall 
intensities, and (2) that the attenuation of all the lossy media (except 
oxygen and water vapor) vary proportionately with the attenuation through 
rainfall. The first of these assumptions may be quite representative of 
the true state of affairs for precipitation from stratiform clouds, except 
that the melting level is lower in the winter months or in cooler climna- 
tological regimes, The second approximation rather crudely attempts t 
account for the variation of cloud and snow content with varying rain 
intensity. 

It was noted earlier that a rainfall rate of 4.5 mn/hr is exceeded 
Curing a total of about 10 hours in the course of a year in England. 
Similarly, rainfalls of 2.6 mm/hr and 0.9 mm/hr are exceeded during a 
total of 1 day and 1 week per year, respectively. Accepting the first 
assumption made above, that the model of Case VI is sufficiently rep- 
resentative of all stratiform rain systems, then one may also estimate 


that Case VI losses are exceeded for about 10 hours ver year. The 


attenuations expected with the lower intensities are found by modifying 


*Poth of these papers deal with the problem of attenuation of microwaves 
in the atmosphere. They should be studied by all concerned with the effects 
of weather on radar. 
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the Case VI losses in accordance with the second approximation. ‘In this 
manner we have made a rough estimate of the total time per year that losses 
of a certain magnitude might be expected under rainy conditions. The 
results are shown in Table 2.2. 
Table 2.2 Two-way attenuation (over a 12-mile vertical path) expected 
to be exceeded for 10 hours, 1 day, and 1 week per year under rainy conditions 


in England. 


Rainfall Rate Time Exceeded Two-way Attenuation (db) 
(mm/hr) (per year) 302 cm 1.8 cm 1.25 cm 0.9 cm 
Ao5 10 hours 6.580) feso. 5.86 7603 
axe 1 day 6.385 Ty 68 ere 
0.9 1 week Oech yale CLD 2-56 Pa aly 


It should be noted that the losses associated with a rain of 0.9 mm/hr 
are roughly similar to those obtained in Cases IV and V, for a heavy 
snowstorm and a dense water cloud, respectively. Thus the total time per 
year during which such losses might be expected should include the time per 
year during which conditions represented by Cases IV and V would be expected. 
Since no satisfactory statistics are available for these conditions, we may 
conclude only that the losses in the bottom row of Table 2.2 are exceeded for 


a minimum of 1 week per year. 
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3. GENERAL SUMMARY 


3,1 Comparison. of Cases 


A complete summary of all the data is given in Tables 3.1, 3.2 and in Fig- 
ure 3.1. ‘The general trends are shown better in the illustration; a more detailed 
description of the trends and their deviations is better obtained in table 3.1. 
Except for Case III, the attenuation tends to increase with increasing case 
number, lHxcept for 1.25 cm at the lower case numbers, attenuation increases 
with decreasing wavelength, Taking the ratio of the total loss at each wave= 
length relative to that at 3.2 cm and averaging over all the cases, we find 
that the atmosphere is about 2.5 times as opaque to 1.8 cm, and about 7 times 
more opaque te 1.25 and 0.9 cm, as it is. to 3.2 cm, It is interesting to 
note, however, that 1.25 cm is better (less lossy) than 0.9 cm in vases IV 
to VII. The use of radar under such weather conditions is, of course, more 
critical since optical techniques fail completely, from the viewpoint of at— 
mospheric losses alone, therefore, we conclude that "the shorter the wavelength, 
the less transparent the atmosphere." This is in contradiction to the re aston 
which evidently led to the selection of 0.9 cm as a replacement for 1.25 cm. 
Except in those cases in which water vapor contributes more than roughly 
50% to the 1.25 cm loss, the atmosphere is more opaque to 0.9 em than to 
1.25 cm, | 

The smallest attenuation to be expected at any wavelength is that 
associated with a dry continental polar (cP) air mass, Among the seven 
models used in this study, Case III comes closest to the clear cP air mass, 

This is the case of an ice crystal cloud at the transition zone between a 
mixture of continental and maritime polar air in the lower layers and cP 


air aloft. The presence of the ice crystal cloud of relatively high density 
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The summary of total attenuation (0-12 miles) versus 
wavelength for the various cases. 


10 20 30 40° 50 60 
2 WAY ATTENUATION - db (RANGE 50MI. X =1.8cm) 


Fig. 3.2 
The variation of total attenuation, by cases, to a ground 
target at a horizontal range of 50 miles as a function of radar 
altitude, 
40 


TABLE 3el . Summary of two-way losses over 12 mi vertical path broken down 
by elements. 
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is overbalanced by the decreased eacer content in the cold air. ven the 
dense ice crystal cloud and heavy snowfall of Case IV do not contribute 
greatly to the total loss except at the shortest wavelength. In that case 
(Case IV, 0.9 cm) loss in ice and snow together is just approaching that due 
to oxygene The moderate supercooled water cloud of Case IV is overwhelmingly 
more effective in decreasing the transparency of the atmosphere than even 
heavy snow and, except at 3.2 cm, loss in the water cloud exceeds that due 
to oxygen, The water cloud of Case V is twice as dense as that of Case IV 
but only slightly exceeds it in total effect because of its higher temperature, 
This water cloud equals the importance of the vapor content of the standard 
(mP) air at 1,25 cm and is much more effective than water vapor at the other 
wavelengths, 

Looking at Table 3,1 from another viewpoint, oxygen contributes less 
than about 0.5 db at all wavelengths, while water vapor adds less than 
0,3 db except under the most extreme conditions or at a wavelength of 1,25 
cm, Even a combination of dense ice cloud and heavy snowfall fails to impair 
the situation by more than about 0.2 db except at 0.9 cm where it may be as 
mach as 0.5 db. 

By far the most lossy media are liquid clouds and precipitation, Cold 


(-L0°C) water clouds of moderate to heavy density (about 0.2 graf m? ) and large 


vertical extent (about 3 miles) can be counted upon for the major portion of 

the attenuation at all wavelengths in the absence of rainfall. However, rain 
of any appreciable intensity (in excess of about 2 mn/hr) and vertical extent 
(about 2 miles) generally exceeds the opacity of all but the densest water 


clouds, 
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3.2 Cumilative Losses Versus Height 


Table 3.2 is a summary of the heights.at which the cumulative losses 
upward from the ground reach 67% and 90% of the 12-mile total. In every 
case the 67% height is below 5 miles in altitude so that a rough approxi- 
mation to the total loss may be obtained by neglecting everything above 
that height. As clouds and precipitation are introduced, the greatest 
increases in attenuation are made in the lower layers thereby lowering 
the 67% and 90% altitudes. The 67% height is always below the top of the 
clouds, and when clouds are in liquid form and relatively dense, the 90% 
height earns near their tops. In the presence of appreciable rainfall, 
the 67% altitude occurs toward the top of the rainfall region, 

When oxygen contributes an appreciable part of the total loss (Cases 
I — III), the 90% height occurs around 7-8 miles, except at 1.25 cm, 
Because of the importance of water vapor at 1.25 cm, both the 67% and 90% 
heights are appreciably lowered in Cases I-III, When water vapor no longer 
is of great significance at 1.25 cm (i.e, when its contribution is less 
than about 50%), there is good agreement between the heights at 0.9 and 
1.25 cm. In general, the 67% and 90% altitudes lower with decreasing 
wavelength. except..for a minimum at 1.25 cm, 

3.3 The 


It was noted earlier that we have used oxygen attenuation coefficients 


Effect..of Using the Lower Oxygen Attenuation Coefficient* 


which lie between two extreme estimates, Recently, the lower of these 
estimates has been more favorably accepted as the correct one, The use of 
the lower values reduces the total 12-mile loss computed for oxygen alone 
by the following amounts: 0,053, 0.067, 0.175 and 0.222 db, the values 


increasing toward the shorter wavelengths. The net effect on cumulative 


* See Appendix, 
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90% of the 12 mile total. 
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losses for the individual cases is shown in Table 3.3. Of course the 
greatest percentage changes. occur under those conditions in which oxygen 
is important. The-maximum decrease-in total loss is 35% at 0.9 cm in 
Gase III, In Cases IV-VII, in which water clouds and rainfall are 
significant, the maximum decrease is 15% even at 3.2 cm. In general, 
therefore, the total computed attenuation cannot be reduced markedly by 
use of the lower oxygen coefficients; also, the worse the weather, the 
smaller the fractional change: in total loss, 
3.4 Optimum Flight Altitudes 

Finally, it is of interest to consider the question of optimum flight 
altitude, or altitude at which the effects of atmospheric attenuation are 
minimized, In this regard it is also useful to know how much improvement 
may be gained by a change in altitude, 

In Figure 3.2 we have plotted, as a function of the altitude of the 
radar, total two-way attenuation over paths extending from an airborne 
radar to a ground range of 50 miles, for a wavelength of jeter: The optimum 
altitude, as we might have anticipated, is the highest one in every instance, 
This is particularly marked in those cases with major losses in the lower 
layers. In moving the radar to higher levels, the radar beam cuts across the 
high attenuation (lower altitude) layers at a greater angle of incidence, 
thereby reducing the distance of propagation through them. However, in Per 
cases where oxygen is an important part of the total 12-mile loss, as in I and 
III, very little improvement can be realized by an increase in altitude, 

When there is a substantial amount of rain or evada liquid water in the 
lower levels, as in Case VI, the rate at which the total loss falls off with 
altitude is quite large and an appreciable gain results. Thus, in Case VI an 


increase of altitude from 7 to 12 miles increases the transparency of the 
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Comparison of total losses using (a) Intermediate 0- 
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6a23 
6.18 
0.80 


21-8 


1369 
0.30 


41.8 
41.6 
0.40 


65.2 
6h09 
0.30 


path by 6,3 db. (This considers the small loss due to the increased slant 
range). The least efficient altitude at which to fly (if the radar is directed 
at the ground) is just above the most lossy layers, Cases IV, V and especially 
VI on Figure 3.2 illustrate this feature very well. This becomes of possible 
importance in view of the standard practice of favoring flight altitudes which 
are just above the cloud decks. This practice should be avoided whenever 
radar surveillance of the ground is a prime consideration, 
4. CONCLUSIONS 

The transparency of the atmosphere to short microwaves is extremely 
variable, depending primarily upon the cloud and precipitation content of 
the lower layers and to a lesser extent (except at 1.25 cm) on the water 
vapor content of the air mass, Because of the complexity of the vertical 
distribution of the several attemating elements (rain, snow, clouds, water 
vapor, and oxygen) and because of the inadequacy of the statistics on the 
frequency with which they occur in various combinations, no satisfactory 
estimate can presently be made of the amount of time during which losses of 
certain magnitude can be expected, Even the very crudest of approximations 
which employs surface rainfall statistics and a physical model of the vertical 
distribution of the associated clouds, snow and water vapor, yields only a 
rough estimate of the minimum fraction of time during which various losses 
may be anticipated at a particular location. The results of such approxima- 
tion made for England are shown in Table 2,2, 

Two-way losses through the vertical extent of the atmosphere (up to 12 
miles) have, however, been computed for a number of physical models of cloud 
and precipitation systems and for clear maritime polar and tropical air 


masses, The decibel losses for these cases are summarized in Table 3.1 and 


Da 


Figure 3,1. ‘The reader may select one or more of these models and cor- 
responding losses as typical of the average or extreme conditions under 

which his radar may be required to operate, Such an approach may well result 
in a better estimate of typical losses than a statistical survey which is 
based on the necessarily crude use of scanty data. 

The cumulative losses over a vertical path from the ground up to an 
altitude less than 12 miles may be obtained directly from the curves for the 
individual cases (Figures 2.3 = 2.15). The heights at which the 12-mile 
losses of Table 3.1 are reduced by 10% and 33% are shown in table 3.2, Fer 
non-vertical paths through a horizontally stratified and homogeneous at-= 
mosphere, Figure 2.5 gives factors by which the vertical path losses must 
be magnified, 

After consideration of the behavior of microwave attenuation as a 
function of wavelength and of the weather composition of the tronosphere, we 
have attempted to draw some general conclusions, these follow: 

1. in general, the opacity of the atmosphere decreases with wavelength. 
Using an artificial average of the seven model weather situations oi the report, 
we find that the atmosphere is about 2.5 times as opaque to 1.8 cm waves, and 
about 7 times more opaque to both 1.25 and 0.9 em, as it is to 3.2 cm radiation. 
(These figures make no account for the frequency of occurrence of the various 
weather models) ‘a 

2. The only advantage which 0.9 cm has over 1,25 cm radiation is its 
smaller absorption by water vapor, Thus, it is superior to 1.25 cm only in 
those cases involving clear, warm, moist air or light clouds and precipita- 
tion, lkadiation at 0.9 cm is also somewhat more affected by scope weather 


clutter so that even in light rain or snow, ground targets may be obscured, 
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However, the overall advantage of either one of these wavelengths over the 
other is not striking. ‘In cold continental or arctic climates where water 
vapor content is low and water clouds and precipitation cannot be signifi- 
cant, the losses at 1.25 and 0.9 cm are likely to be very much the same, 

3, If one had to select arbitrarily an artificial world-wide average 
of the bad weather conditions in which the use of radar will be required, we 
would guess at something between Cases IV and V (Table 3.1). The two-way 
losses over a vertical path through the lower el 2 miles of the atmosphere 
would then be of the order of 0.4, 0.8, 2,0 and 2.5 db at 3.2, 1.8, 1.25 
and 0.9 cm respectively. Such losses are not unlike those to be expected in 
a light rain having a vertical structure similar to that of Case VI (spe 
Table 2.2). Over an oblique path from an altitude of 8 miles to a ground 
range of 50 miles, these values would be magnified by a factor of about 6.6 

4. Clouds and precipitation in the liquid state are the most important 
attenuating media in the troposphere, except that oxygen is relatively of 
substantial importance at 3.2 cm, and water vapor at 1.25 em, Cold (~-10°C) 
water clouds of moderate to heavy density (about 0.2 gm/m? ) and large 
vertical extent account for the major portion of the losses at all wave=- 
lengths in the absence of rainfall, However, rain of any appreciable inten- 
sity (in excess of about 2 mm/hr) and vertical depth generally exceeds 
the opacity of all but the densest water clouds. 

52 Two-thirds of the total 12-mile loss generally occurs in the lower 
5 miles of the atmosphere or, when clouds are present, below the cloud tops. 

6, In using airborne radar for ground surveillance, the worst leveliat 
which to fly is just above the most lossy layers (generally just above the 
cloud tops) . Substantial gains in oblique transmission may be had by in- 
creasing flight altitude when the lower layers contribute an appreciable loss, 


The practice of flying just above cloud decks should therefore be avoided. 
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APPENDIX & 


Oxygen Absorption 

The variation of oxygen absorption with altitude which we have used 
throughout Part II of this study is based upon a statement by Van Vleck 
(1951, pe 647) to the effect that the oxygen véetitoiene varies as the 
quotient of density by temperature, or directly as the pressure and inversely 
as the temperature squared. However, this is evidently incorrect, The 
proper variation with altitude should include the square of the pressure as 
a factor because the wavelengths of interest fall on the wing of an absorption 
band (see Section 7.2, Part I and Kogers, 1951). Thus, the oxygen curves of 
Figure 2.2 (amd the remaining illustrations) should fall off more rapidly with 
altitude than shown, Rogers (1953) has computed the cumulative losses by 
oxygen and water vapor as a function of altitude using the proper pressure 
variation and the highest estimate of the oxygen absorption seat eieiennee At 
an altitude of 7 miles, his cumulative losses are very nearly equal to those 
which we have computed using the intermediate estimate for oxygen and the 
slower decrease with altitude, Averaged over the lower 7 miles of the atmos— 
phere, therefore, we have effectively used the higher estimate of the oxygen 
coefficient, This is very fortunate and consequently, we have seen no need 
to revise our results al though the oxygen losses in the layer from 7-12 miles 
will be slightly on the high side, 

In summary, although attenuation by oxygen should decrease more rapidly 
with altitude than shown in our figures, the cumulative losses at least up | 
to an altitude of 7 miles are no larger than those to be expected with the 


highest estimate of the oxygen absorption coefficient, 
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APPENDIX B 


The performance of a weather radar may be considered under headings of 
resolution and sensitivity. Decreasing the wavelength would increase both 
these factors if it were not for an accompanying increase in attenuation by 
rain, The attenuation works against the enhanced resolution by introducing 
distortion, nearby precipitation patterns casting shadows on more distant 
patterns. It works powerfully against the enhanced sensitivity; at shorter 
wavelengths sensitivity depends more on the amount of intervening rain than 
it does on the distance away of the target rain, The present report will not 
deal with resolution, but will deal carefully with the matter of sensitivity. 


The performance of a weather radar may be described by the equation 
Po h ARS 
Py = Mae ( Age), (1) 


where P, is the power received at the radar, P, is the power transmitted, 

h is the pulse length, A is the geometrical area of the antenna, 36> is 
the sum of the back-scatter cross-sections of all the precipitation particles 
in unit volume, andr is the range. Ae is the effective area of the antenna 
for extended targets and is less than A. The exact value of Ae is not 

easily determined. For this reason A only is retained in the equation and a 
conservative estimate (1/2) is made for Ae/A. The term F is a factor due to 
unknown causes which a careful experimental check by the M.I.T, Weather Radar 
Project has shown to be about 1/5. Therefore we have allowed a value of 0.10 
for factors in the bracket. It is felt that with these factors, Eq. (1) re- 
presents the performance of a well-maintained radar as closely as present know- 


ledge permits. 


Equation (1) is valid, provided the beam is filled with scatterers and 
provided that attenuation due to intervening rain and atmospheric gases may 
be neglected. 


The back-scatter cross section of a raindrop is given, according to 
Rayleigh's approximation, as 


> 6 2 
& = ee ; (2) 
PN Be ie 1 
m- + 2 
tive index of water, and is the wavelength of the radar, Rayleigh's 
approximation is applicable (see Fig. 5x, Marshall, East and Gunn, 1952). 


where D is the drop diameter, WR > m being the complex refrac- 


The sum of the sixth powers of the drop diameters in unit volume is 
given by 6 6 
ZED? = 2.0 RLeS yQ710 om? om (3) 


where R is the rate of rainfall in mm ho (This relation is also quoted 
from Marshall, et al (1952). ) It depends on the size distribution of the 
drops, as obtained by Laws and Parsons (1943) and Marshall and Palmer (1948). 


Combining Eqs. (1), (2) and (3), and solving for R one obtains: 
: 262 ee +0025 
Be er Pete Leo sere alhlaem ye 


<All 
with all the quantities in c.g.s. units, and R inmmhr . (Here A,/AF has 
been put equal to 0.10 as discussed above.) 7 


The minimum detectable rainfall R, may be calculated from this equation 
if, for Pp, one uses the minimum detectable signal power, At the 1952 Radar 
Weather Conference such minimum values were reviewed (see Appendix Tt); and the 
consensus of opinion seemed to be that min. detectable Pr is about 4 x 10713 
watts for both 3.2 and 5.7 em radar. This is the value of Py that will be used 
throughout this report. 


Putting 1,4 be = 0.93 (a value very nearly independent of wavelength) we 
rewrite Eq. (4) for the minimum detectable rainfall in more convenient units: 


2eD0,1.25 P 0.625 P : —0,625 —125 
Pa d. us PE eee 4) 
Rm = 0.0981 ( 57 =) ( hx10-13 es fs kw 300 7 _ 


(1am es (mi) (5) 


where d is the diameter of the (circular) antenna, 


When the beam is not filled by pecipitation at range r, the equation 
needs to be modified. The diameter of a circular beam is (between hal f- 
power points) r 1.14 ,/d; hence, if we view a storm of linear dimension iy. 
it follows that the beam will be just filled at range rg, such that 


Po eee A= L, (6) 


Beyond ro the beam will no longer be filled, and here P, in Eq. (1) should be 
proportional to r4 rather than to r-<, Carrying through this modification, 
taking care that the two forms of Eq. (1) mst agree at r = Yo, results in 


; seid a 50 2 Pal P ) 625 
a A re r fot Parise hae : OR 
Ri = 0.00218 ( = =) ( epee -) Pa cate 306 n 


(mm hr~) (mi) = -1, : 
ee | es 


Equations (5) and (7) give the minimum detectable rainfall without allow- 
ing for attenuation by intervening rain and by air and water vapor. Such 
attenuation causes a drop in the received power, and may be allowed for most 
easily by a corresponding drop in P,, or an equivalent increase in Rie Mees 
the above equations. 


Attenuation by rain may be expressed accurately by 


radar attenuation by rain = oe fa R® dr, (8) 
(in db) - @) 


where 


202 
0.0288 


1.3 


Table I 


provided R is in mm hr 1 and r in miles (data from Marshall, et al (1952) ). 
The values of & are sufficiently close to unity to allow the simplification: 


radar attenuation by rain = K x quantity of Sree rain 
(in db) in mm hr- x miles) (9) 


where K is only slightly dependent on R. If R is about 30 mm hr=l, 
K has the values: 


= ‘ea Coe oe Pa 


0.013 db/(mm hr-l mi) 


Table 2 


3 


The attenuation by atmospheric gases (notably oxygen and water vapor) 
varies with temperature, pressure and humidity. Since the values are small, 
only representative figures, corresponding to a temperature just above 9°C, 
standard pressure, ani water vapor pressure close to saturation, need be. 
used, Thus, the radar attenuation by gases, k, may be taken to be 


E | 3.2 | 
O 


| 5.7 cm 
r 


k | 2054 | 0,036 db ni | 


: 


Table 3 
These values are taken largely from Ryde (1946). 


There are, of course, a large number of ways in which the variation of 
the minimum detectable rainfall with the various radar parameters, with in- 
tervening attenuation and with the extent and intensity of the storm to be 
detected, may be displayed. One practical and useful way is to show the 
intervening rain y (in mm hr7+ x miles) , through which a storm of given ex- 
tent and intensity may be detected, as a function of range, r. Neglecting 
attenuation by the gases, this relationship may be derived as follows. (It is 
sketched in dashed lines in Fig. 1) 


y 


=4 


interven ing vain Cmm hr mi) 


/000 ———— without allowing for 


gas attenuation 


eee lowing for gas 
attenuation 


500 


§0 
range r (statute miles) 


Au eel 
aby 


The radar detects Pals) (the minimum detectable rainfall at range r, 
as given by Eqs. (5) or (7)). If one is only interested in detecting 
Ro (a more intense rainfall), then the difference in power 


leds 
10 log te Be 


may be used up by attenuation, Thus 


166 
Ky = 10 log ( E2 ° (10) 


As long as the beam is filled with precipitation, Eq. (5) applies for 
Rm, and 


20 16 . 
y=-f logr+ K 108 Cat (11) 


Beyond range r,, the beam is no longer filled and instead of Eq. (11) one 
has (combining Eqs. (10) and (7): 


16 R 
yin WO logr+ > logge (12) 
yess he Rim . 
Here as miges and © 72650 » (13) 


are constants depending on the radar parameters, which may be obtained from 
Eqs. (5) and (7). Plots of y and y' against log r are straight lines. 


Allowance for gas attenuation is easily made, simply by reduc in valor 
value of r, the value of y (or y') by an amount 


Ay= Ayt= 2, (14) 


This correction is quite small, so that in spite of it, the plots of y or 
y' vs. log r still appear to be nearly straight (the solid lines of Fig, 1 
and all lines of Fig. 2). 


The effect of wavelength is emphasized here, Loci of y or y' are 
therefore drawn in pairs: one locus for 3.2 and one for 5.7 cm equipment, 
with all other parameters unchanged, Atypical example of such a pair 
is sketched roughly in Fig, l. 


Of special interest in the point of intersection (I) of the two curves. 
The ordinate (y,) of this point is the intervening rain at which the two 
wavelengths have the same range. For lower values of y, the shorter wave- 
length has the better range; for higher values of y, the longer wavelength 
is much to be preferred. Neglecting gas attenuation, y, turns out to be 
2ndependent of all radar parameters and of Rp and L, provided the beams are 
either both filled or both unfilled. If the 3.2 cm beam is filled, but the 
5,7 cm beam is not (a frequent case), y. turns out to depend somewhat on trans- 
mitter power, pulse length, etc., but not on the antenna size. The effect on 
Yo of the attenuation by the gases is generally unimportant. 


Figure 2 shows the final results of the calculations in the form out- 
lined above; construction lines have been omitted and three pairs of loci are 
given, each pair having. one member of 3.2 cm and one for 5.7 cm. This chart 
may best be discussed in terms of a particular pair of loci: the solid lines, 
for instance, refer to transmitter power 120 kw, antenna aperture 18 inches 
in diameter, pulse length 300 meters, intensity of target rainfall 10 mm/hr. 

For amounts of intervening rain up to 200 mm hr-1 miles, the 3.2-cm equipment 
gives greater ranges than the 5.7 cm. With this much intervening rain, both 
wavelengths have a limiting range of 31 miles, As one proceeds to greater 
amounts of intervening rain, the range at 3.2 cm continues to drop off very 
rapidly, so that with twice this much intervening rain the range is about 

7 miles, At 5.7 em, on the other hand, the same 400 mom/ hr of rain drops the 
range back by less than 15 per cent to 27 miles, It was necessary to specify 
the extent of the target shower, since this determines whether or not the target 
precipitation completely fills the cross section of the radar beam, A shower 

of linear dimensions 3 miles high by 3 miles wide was chosen rather arbitrarily . 
for all curves, At 5.7 cm this target fills the beam of the 18" antenna at 
ranges up to 21 miles, the narrower 3,2 cm beam at ranges to 38 miles, Kinks 

in the loci under discussion may be noted at those ranges. 


The amount of intervening rain at which the two wavelengths have the 
same range, i.e., about 200 mm hewt mi, is not very much in stormy weather 
situations, and so attenuation at 3.2 em is seen to be a very serious diffi- 
culty. It is not just a matter of the range being limited, there is the un- 
certainty whether one sees light rain through a small amount of intervening 
precipitation or heavy target rain through much intervening precipitation. 


Let us now look at the altermatives. Reduction of the power by a fac— 
tor 3 or of the minimum detectable target rainfall to 5 mm brvL makes for a 
serious further reduction in range (dotted lines). On the other hand, the - 
effect of increasing the size of dish is to extend the range very consider- 
ably (broken lines). The range at the longer wavelength comes between 
57 miles in the clear, ami 45 miles with about 400 intervening mm hr= miles. 
Effectively the same result can be achieved by increasing the target rain-—. 
fall to 50 mm hr7+, | 


The chart shows three pairs of curves covering four sets of parameters 
at the two wavelengths. It is, however, easy to construct further loci by 
moving the lines, without changing their direction, in the following way. 
(Such an extension of the chart is only approximate, as this procedure should 
really be applied to the construction lines of Fig. 1.) 


To allow for variation in 


(a) transmitter power P,, move curves vertically, changing the height of the 
kink (or any other point of fixed range) by 


yt -~y-% x (change in Pp in db) 


1 
(= = 10g Po'/Po), 


(b) pulse length h, proceed exactly as for variation in P,, 


(c) target rainfall Ro, move curves vertically, changing the height of 
the kink, etc., by 


yi-y= iy (change in Ry in db) 


(= = log Ry /Ro) 


(d) antenna diameter d, move curves horizontally, adjusting the range 
of the kink proportionately to d, (ise. ro'/ro = d'/d); 


(e) linear dimension of target rain, L, move kink along left part of 
locus (produced, if necessary) to a range proportional to L; (i.e. Se ./te et 
L!/L). Then draw right part of locus parallel to the right part of the 
locus initially given. 


It is coneluded from this investigation that attenuation by intervening 
rainfall effects a very serious limitation to the sensitivity of 3 cm weather 
radar. A frontal line of showers is likely to involve the passage of the 
radar beam through a few hundred mm hr7! miles of rain, probably, though not 
certainly, less than 500. Sensitivity at 3 cm varies very rapidly with the 
amount of intervening rain, and this rapid variation in itself would lead 
to uncertainty in interpreting signals. Five hundred mm hr7- miles is enough 
to render 3 cm equipment practically inoperative, 


At wavelength 5.7 cm, attenuation is appreciable but the sensitivity 
changes more gradually with intervening rain, and the range would never be re- 
duced by more than 30 per cent. (This assumes that the shower does not fill 
the beam; in all likelihood this will be the case at limiting ranges with the 
relativeity small size of an airborne antenna.) At this wavelength, however, 
it will be difficult. to achieve the maximum desired ranges even in clear air. 


The most promising solution would appear to be the use of wavelength 5.7 em, 
doing everything possible in design and construction to keep up the sensitivity. 
At the same time, it may be reassuring to note that the range increases with 
the, intensity of target rainfall. The range is doubled in going from the 10 mm 
hr specified here to 60 mm hr-l, Radar observations on thundershowers almost 
always indicate intensities of this order at the core, or at any rabe a core 
condition providing a signal equivalent to 60 mm hr-l in radar reflectivity. 


The two wavelengths have been compared only in sensitivity. For a 
constant size of antenna, they will also differ in resolution: the shorter 
wavelength will have a correspondingly narrower beam and correspondingly 
higher resolution. But attenuation will introduce distortion into the pic-— 
tures, altering outlines and casting shadows of nearer showers on more distant 
ones. t may be that the initially better resolution of the shorter wave- 
length will be more than cancelled out by attenuation, in the same way as the 
initially higher sensitivity. The situation with regard to resolution and 
distortion has not been studied; the best reference to date woullappear to be 
Atlas and Banks (1951). : 


TABLE 4 


Performance of Radars of Different Wavelengths 


This table is taken from a brief note (pages B-37/38) in the Proceedings of 
the Third Radar Weather Conference, held at McGill University, 17-19 September 
1952. The data represent the general consensus of the Conference, which com— 
prised the majority of radar weather specialists. 


N (cm) 0.9 1.25 3,2 5.6 10,0 


Maximum available 20 1,0 350 300 1000 
transmitter power (P,,kw) 

Minimum receivable wo 5 x108 yxio’ 3x10B 2x10 
signal power (P,, watts) 

Performance factor 

(P,/P,x ,-*, relative ( 61 65 17 20h 1 

to 10 cm equipment) 

Performance level (db) 18 18 12 3 ) 


Attenuation (two-way transmission) 


Cloud (density 1 gm n)t ae 1.75 0,28 

(db mi-l) 

Rain - db (mm hr™) x mi)74 0,09 0.015 0,00 
(at R = 100 mm hr@+) 
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